Classification and Properties of Symmetry Enriched Topological Phases: A 
Chern-Simons approach with applications to Z2 spin liquids 
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We study 2+1 dimensional phases with topological order, such as fractional quantum Hall states 
and gapped spin liquids, in the presence of global symmetries. Phases that share the same topological 
order can then differ depending on the action of symmetry, leading to symmetry enriched topological 
(SET) phases. Ilere we present a K-matrix Chern-Simons approach to identify all distinct phases 
with Abelian topological order, in the presence of unitary or anti-unitary global symmetries . A 
key step is the identification of an edge sewing condition that is used to check if two putative 
phases are indeed distinct. We illustrate this method for the case of Z2 topological order (Z2 spin 
liquids), in the presence of an internal Z2 global symmetry. We find 6 distinct phases. The well 
known quantum number fractionalization patterns account for half of these states. Phases also 
differ due to the addition of a symmetry protected topological (SPT) phase. Also, we allow for the 
unconventional possibility that anyons are exchanged by the symmetry. This leads to 4 additional 
phases with symmetry protected Majorana edge modes. Other routes to realizing protected edge 
states in SET phases are identified. Symmetry enriched Laughlin states and double semion theories 
are also discussed. Two surprising lessons that emerge are: (i) gauging the global symmetry of 
distinct SET phases can lead to states with the same topological order (ii) gauge theories with 
distinct Dijkgraaf-Witten topological terms may have the same topological order. 
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I. INTRODUCTION 

It was long believed that phases of matter arose from 
different patterns of symmetry hieakin^^. The discov- 
ery of integer"^ and fractional^ quantum Hall (FQH) ef- 
fects demonstrated however that there exist many differ- 
ent phases of matter which lie outside this paradigm. In 
particular, the FQH states differ in their internal quasi- 
particle structure as we ll as their boundary excitations 
and are said to realizd^ 'topological order' which implies 
ground state degeneracy (GSD) on a closed manifold (a 
Riemann surface of genus g) and emergent anyon excita- 
tions which obey neither bosonic nor fermionic statistics. 
Another class of topologically ordered phases are gapped 
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quantum spin liquid^. Recently, several examples of 
gapped spin liquids have appeared in numerical calcula- 
tions of fairly natural spin 1/2 Heisenberg models, on the 
Kagome^and square lattice (with nearest and next neigh- 
bor exchange)^ ^. Calculation s of e ntanglement entropy 
point to Z2 topological ordei'^^'^. However, the pre- 
cise identification of these phases require understanding 
the interplay between topological order and symmetry 
in these systems. The symmetries include both on-site 
global spin rotation and time reversal symmetries, as well 
as the space group symmetries of the lattice. Kagome lat- 
tice antiferromagnets, such as herbertsmithite, may pro- 
vide experimental realization of this physics, although ex- 
perimental challenges arising from disorder and residual 
interactions continue to be actively studied. This moti- 
vates the study of distinct topologically ordered phases 
that may arise in the presence of symmetrjE^tl^ 

In the presence of symmetry, the structure of topo- 
logical orders is even richer. The microscopic degrees of 
freedom in the system are either bosons or fermions, and 
they must form a linear representation of the symme- 
try group Gg. The emergent anyons, however, doesn't 
need to form a linear representation of Gs- Instead they 
could transform projectively under symmetry operation, 
i.e. each of them can carry a fractional quantum number 
of symmetry. For example, the Laughlin FQH states at 
filling fraction v = \/rn^, where each elementary quasi- 
particle carries a fractional (l/m) of the electron charge. 
This phenomena is widely known as fractionalization, al- 
though a more approp riate name is perhaps symmetry 
fractionalizatioT^^^^^. And the associated symmetry in 
Laughlin states is the C/(l) charge conservation of elec- 
trons. While the emergent quasiparticles transform pro- 
jectively (instead of hnearly) , the microscopic degrees of 
freedom always transform linearly under symmetry, sim- 
ply because each microscopic degrees of freedom can be 
regarded as a conglomerate of multiple emergent quasi- 
particles. 

Even in the absence of topological order, when symme- 
try Gs is preserved, different symmetry protected topo- 
logical (SPT) phased^ emerge which are separated from 
each other through phase transitions. These SPT phases 
feature symmetry protected boundary states which will 
be gapless, unless symmetry Gs is (spontaneously or ex- 
plicitly) broken on the boundary. Well kno wn ex am- 
ples of SPT phases are topological insulatord^HHl and 
superconductors^'^. In 2+1-D all SPT phases have sym- 
metry protected non-chiral edge modeP^^t^Sl, 

The existence of SPT phases further enrich the struc- 
ture of topological orders in the presence of symmetry. 
In other words, topologically ordered phase is not fully 
determined by how its (anyon) quasiparticles transform 
(projectively or not) under symmetry: its microscopic de- 
grees of freedom could form a SPT state in parallel with 
the topological ovde^^. The formation of SPT state will 
e.g. bring in new structures to the edge states of the 
topologically ordered system, and lead to a distinct sym- 
metry enriched topological (SET) order. Therefore, two 
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Only microscopic degrees of freedom ("electrons") 
can tunnel between the two edges! 



Figure 1: (color online) Edge Sewing Criterion to distinguish 
symmetry enriched topological (SET) phases. Only the mi- 
croscopic degrees of freedom i.e. "electrons" (and not gauge 
charged objects such as anyons/fractionalized quasiparticles) 
can tunnel between the two edges of a pair of semi-infinite 
cylinders. If two SET phases can be continuously tuned into 
one another without a phase transition (while preserving sym- 
metry), there is a "smooth" sewing between the two cylinders 
of SET phases #1 and #2. This implies that all edge ex- 
citations are gapped by a few symmetry-allowed terms that 
tunnel "electrons" between the two edges. In the thermo- 
dynamic limit these tunneling terms lead to M degenerate 
ground states, corresponding exactly to the M-fold torus de- 
generacy of the topological order. On the other hand, if 
the two SET phases are different, there is no such "smooth" 
boundary condition to sew the two edges. A precise version 
of this statement is formulated in Criterion I in Section fll El 



different SET phases sharing the same topological order 
can differ by the symmetry transformation on their anyon 
quasiparticles, or by their distinct boundary excitations. 
Clearly we have a question here: given two states shar- 
ing the same topological order while preserving symme- 
try Gs, can they be continuously connected to each other 
without a phase transition, if symmetry Gs is preserved? 

We address this issue for 2+1-D (Abelian) topological 
orders. Focusing on on-site (instead of spatial) symme- 
tries, we present a universal criterion (Criterion I in sec- 
tion II El related to the edge states of these 2+1-D SET 
phases, which works for both unitary and anti-unitary 
on-site symmetries. The physical picture behind this cri- 
terion is demonstrated in FIG. [l] Two SET phases #1 
and #2 living on the two cylinders are considered the 
same if they can be smoothly connected together via tun- 
neling of microscopic degrees of freedom between the two 
edges. Distinct SET phases on the other hand, present 
an obstruction to such a smooth sewing. 

The above criterion allows us to clarify the struc- 
ture of symmetry enriched topological (SET) orders in 
2+1-D. The method we follow is the Chern-Simons ap- 
proach, which provides a unified description for low- 
energy bulk and edge properties of a generic Abelian 
topological order^3i 2ll in 2+1-D. In particular the bulk- 
edge correspondencies'^ in Chern-Simons approach en- 
ables us identify all edge excitations with their bulk 
counterparts, such as the microscopic degrees of freedom 
(bosons/fermions) and anyons. Therefore the above cri- 
terion of smooth sewing boundary conditions for two dif- 
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ferent SET phases can be made precise within the Chern- 



Simons approach (see section HE). 

More concretely, a 2+1-D Abelian topological phase is 
fully characterized by a symmetric integer matrix K in 
the Chern-Simons approach. When symmetry Gs is pre- 
served in the system, the anyons could carry a fractional 
symmetry quantum number (or transform projectively 
under the symmetry^^), while the microscopic degrees of 
freedom (bosons/fermions) must form linear representa- 
tions of the symmetry group Gg — {g}. The relation be- 
tween microscopic degrees of freedom and fractionalized 
anyon excitations is especially clear in the Chern-Simons 
approach. 

Based on the mentioned Criterion I to differenti- 
ate distinct SET phases, we can classify all different 
SET phases with the same topological order {K} 
and symmetry G^. In this work, we studied va rious 
examples: Z2 spin liquids^^, double semion theorjl^^lMI 
and bosonic/fermionic Laughlin states- at filling fraction 
1/ — 1/m. We consider both anti-unitary time reversal 
symmetry Gg = ■, and unitary Gg = Z2OVZ2 x Z2. 
These are the analogs of the spin rotation symmetry of 
Heisenberg magnets. Classification of these SET phases 
with symmetry Gg are summarized in TABLE |l]|V] 



Unconventional SET Phases: We divide SET phases 
into two types - conventional and unconventional. In con- 
ventional SET phases, all (anyon) quasiparticles merely 
obtain a J7(l) Berry phase under any symmetry oper- 
ation. In contrast, in the more exotic 'unconventional' 
SET phases, certain symmetry operations exchange two 
inequivalent anyons, instead of just acquiring C/(l) Berry 
phases. For example, under on-site unitary Z^ symmetry 
operation, the two anyons i.e. the electric charge e and 
magnetic charge m of a spin liquids are exchanged (see 
TABLE III). Previously, such a transformation law was 



double semion models, do not have gapless excitations 
at the edge. However, these may appear with additional 
symmetry. Indeed, the 'unconventional' Z^ SET phases 
have Majorana edge states. Since they are protected 
by an onsite Z2 symmetry, they are stable even in the 
presence of disorder that breaks translation symmetry 
along the edge. Two further mechanisms for gapless edge 
modes in 'conventional' SET phases may be identified. 
The first is the trivial observation that adding an SPT 
phase leads to a corresponding protected edge state. The 
second mechanism operates when both the electric and 
magnetic particle of the Z2 gauge theory transforms pro- 
jectively under symmetry. Then, one cannot condense 
neither of them at the edge - implying a protected edge. 
Details appear in Sec |IIIF[ 

Gauging Symmetry: A powerful tool in studying the 
effect of an onsite unitary symmetry G g is the conse- 
quence of gauging i t^^ l ^^ l This means the global Gg sym- 
metry is promoted to a local "gauge symmetry", which 
leads to new topological orders. Distinct topological or- 
ders can help distinguish different actions of the sym- 
metry in the ungauged theory. By this procedure in 
2+1-D, nonlinear sigma models with topological terms, 
which describe SPT phaseJ^SI c an be mapped to gauge 
theories with a topological terrrP^H] Jigcussed by Di- 
jkgraaf and WitterPsl In this work we systematically 
study the consequences of gauging unitary on-site sym- 
metry in Abelian SET phases. It turns out for "conven- 
tional" (type-I) SET phases, the new topological order 
obtained by gauging symmetry is always Abelian and 
Chern-Simons theory is a natural framework to derive 
it. In contrast to the belief that different 'conventional' 
SET phases always lead to different topological orders 
after gauging its symmetry, we found that many differ- 
ent SET phases (with a common topological order and 
symmetry group Gs) can lead to the same topological or- 
der once the symmetry is gauge cPl, as shown in TABLE 



considered in the Wen 'plaquette' modepSH221 for trans- 
lation symmetry, in contrast to the internal Z2 symme- 
try considered here. These "unconventional" SET phases 
have some striking properties. First, the edge features 
gapless Majorana edge modes that are protected by sym- 
metry. Next, if Z2 symmetry is broken at the edge, then 
a Majorana fermion is trapped at the edge domain wall. 
Finally, as illustrated in FIG. |2] when a pair of electric 
charge (e) is created at opposite sides of a sphere, we can 
divide the system into two subsystems A and i?, so that 
there is one electric charge e localized in each subsys- 
tem. Now if we perform the Z2 symmetry operation only 
in subsystem A (fiip all the spins), the electric charge 
e therein will become a magnetic charge to. Since an 
electric charge e and a magnetic charge to differs by a 
fermion f (e x f — m or m x f — e) in the Z2 spin liq- 
uid, this means a fermion mode / must simultaneously 
appear at the boundary separating subsystem A and B, 
as the Ising symmetry is acted on A. This is discussed 
in Section lllT¥2l and Section UTTbH 

Symmetry Protected Edge States: In general, the non- 
chiral topological orders, like Z2 topological order and can be introduced to extract all information of the non- 



[n] and IV Somewhat surprisingly, this furnishes exam- 
ples where a gauge theory with two distinct Dijkgraaf- 
Witten'^^ topological terms, correspond to the same topo- 
logical order. Here the topological terms arising for the 
gauge group Z2 x Z2 are obtained by gauging SPT phases 
and correspond to elements of H^{Z2 x Z2, U{1)). The- 
ories for distinct elements are shown to be equivalent on 
relabeling quasiparticles (an SL{A,Z) transformation). 
Therefore the distinction between these theories requires 
additional information such as specification of electric vs. 
magnetic charges (Section [HI B 3 1 . 

For "unconventional" (type-II) SET phases, however, 
gauging the symmetry leads to non-Abelian topolog- 
ical orders. For example the unconventional Ising- 
symmetry-enriched Z2 spin liquids, after gauging the 
Ising {Gg = Z2) symmetry, lead to non-Abelian topo- 
logical orders with 9-fold GSD on a torus. Interestingly, 
they can be naturally embedded within Kitaev's 16-fold 
way clas sific atio rP^ of 2+1-D Z2 gauge theories (see TA- 
VI I . In this case a vertex algebra approach^ 
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Abelian topological order. In particular after gauging 
the on-site Ising {Z2) symmetry, a new quasiparticle qg 
(called Z2 vortex) emerges as deconfined excitations. It 
is a non-Abelian anyon in the unconventional (type-II) 
SET case, which corresponds to the edge domain wall 
bound state in FIG. HI 

Spin 1/2 From K-Matrix CS Theory: We demonstrate 
how an emergent 'spin 1/2' excitation can be realized 
in the Chern Simons formalism, by studying Z2 gauge 
theories with Z2 x Z2 symmetry. The latter has a projec- 
tive representation that can protect a two fold degenerate 
state, analogous to spin 1/2. This is accomplished by ex- 
panding the 2x2 K-matrix of a Z2 gauge theory to a 4 x 4 
matrix by adding a trivial insulator. Symmetry transfor- 
mations implemented in this expanded space have the 
desired properties (in Sectior III E I . 

Connection to Other Work: A symmetry based ap- 
proach was used to classify Z2 spin liquids in RefflSl 
An advantage of that approach is that it treated 
both internal and space group symmetries. However, 
topological distinctions and the appearance of edge 
states are not captured. Also, the 'unconventional' 
symmetry realizations were not discussed. Finally, as 
mentioned in Ref, 15 the symmetry based approach 
produces forbidden Z2 SETs, that cannot be realized in 
2D, but only as the surface state of a 3D SPT phas^^l 
Our K-matrix approach does not produce such states. 
A different classification scheme in Ref 16, produces a 
subset of our 'conventional' phases although explicit 
lattice realizations are given for them. Finally, RefilTi 
gave a classification based on gauging the symmetry, 
which misses distinctions between phases as discussed 
previously. Our approach is perhaps closest to that 
adopted in ReflHl which however was restricted to 
time reversal symmetric topological states. Thus the 
results in this paper go beyond previous classifications 
of Z2-symmetry-enriched Z2 gauge theories (including 
Z2 spin liquid and double semion theory) , and a detailed 
comparison is given in Section [ill B 3| 



This paper is organized as follows. In Section [IT] 
we introduce the Chern-Simons K-matrix approach to 
(Abelian) symmetry enriched topological (SET) phases 
in 2+1-D. Rules for implementing on-site symmetry in 
a topologically ordered phase are discussed in Section 
|IID| with criteria to differentiate distinct SET phases 
in Section |IIE| Next, in Section |III| we demonstrate 
our approach by classifying SET phases in a few exam- 
ples. They include:(i) Z2 spin liquid with (anti-unitary) 
time reversal symmetry {Gs — Z2) symmetry (section 
mA) TABLE |l|,(ii) Z2 spin liquid with unitary Ising 
(Gr = Z2) symmetry (section |III B[ TABLE In] and lllT| 



(iii) double semion theory with unitary Ising (Gs = Z2) 
symmetry (Appendix [C) TABLE |lV| and (iv) even- 
denominator bosonic Laughlin state with unitary Ising 
{Gs = Z2) symmetry (section |III C[ TABLE |v|). Ap- 
pendix [B] and [D] provide a detailed instruction on how 
to gauge unitary symmetries in 2+1-D SET phase. 



II. CHERN-SIMONS APPROACH TO 
SYMMETRY ENRICHED ABELIAN 
TOPOLOGICAL ORDERS IN 2+1-D 

A. Chern-Simons theory description of 2+1-D 
Abehan topological orders 

In two spatial dimensions, a generic gapped phase 
of matter is believed to be described by a low-energy 
effectiv e Cher n-Simons theory in the long-wavelength 
limiijllHlIIllIlI] Both the bulk anyon excitations and 
the gariess edge states are captured by the effective 

theory^. Examples include integer and fractiona l quan- 
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tum Hall states^, gapped quantum spin liquidg^SHffiil 
and topological insulators/superconductors. When we 
restrict ourselves to the case of gapped Abelian phases 
where all the elementary excitations in the bulk obey 
Abelian statistic^^, a complete description is given in 
terms of Abehan U{1)^ Chern-Simons theory ^'^ ^^ To 
be specific, the low-energy effective Lagrangian of J7(l)^ 
Chern-Simons theory has the following generic form 



N 



i,j=i 



N 
7=1 



(1) 



where fJ-,v, X = 0, 1, 2 in 2+1-D and summation over re- 
peated indices are always assumed. Here • • • represents 
higher-order terms, such as Maxwell terms ^ (e^iyS^af,)^. 
K is a symmetric N x N matrix with integer entries. No- 
tice that the U{1) gauge fields are all compact in the 
sense that they are coupled to quantized gauge charges 
with currents jj. In the first quantized language the 
quantized quasiparticle currents jj are written as 

V/ = l,... ,iV: j?(r)=^4")<5(r-r(")), 

n 

jT(r)=^4"^ri"'<5(r-r(")), a = 1,2. 



(rj^" , r2 ) denotes the position of the n- 



where r^^"^ 

th quasiparticle, and gauge charges Zj"^ are all quantized 
as integers. We can simply label the n-th quasiparticle 
by its gauge charge vector 1^"^ = ,zj^'')-^. The 

self(exchange) statistics of a quasiparticle 1 is given by 
its statistical angle 



ttFk-^, 1 e Z 



N 



(2) 



while the mutual(braiding) statistics of a quasiparticle 1 
and 1' is characterized by 



N 



(3) 



The above statistics comes from the nonlocal Hopf 
LagrangiarP^ of currents jj, obtained by integrating out 
the gauge fields in |r|) . A simple observation from ([s]) 
is that for a quasiparticle excitation with gauge charge 



1 = Kl, 1 e Z 



N 



(4) 
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its mutual statistical with any other quasiparticle 1' is a 
multiple of 27r. In other words, the quasiparticles 1 = Kl 
are locaP^ with respect to any other quasiparticles 1'. 
Therefore they are interpreted as the "gauge-invariant" 
microscopic degrees of freedom in the physical system: 
such as electronics! ^ fractional quantum Hall state, 
and spin-1 magnons in a spin- 1/2 Z2 spin liquicP^. An- 
other direct observation is that when all diagonal ele- 
ments of matrix K are even integers, the microscopic de- 
grees of freedom have bosonic statistics 9 = mod 27r, 
and ([ij describes a bosonic system. When at least one di- 
agonal elements of K are odd integers, there are fermionic 
microscopic degrees of freedom in the system. 

The ground state degeneracy (GSD), as an important 
character for the topologically ordered phase described 
by effective theory ^ is^^ 

GSD = |detK|f. 

on a Riemann surface of genus g. On the torus with 
g = 1, the corresponding GSD= |detK| also equals 
the numbers of differ ent a nyon types (or the number of 
superselection sectorP-Sl—^ in the 2+1-D topological or- 
dered system. A simple picture is the following: any two 
anyons differing by a (local) microscopic excitations are 
the same (or more precisely, belong to the same super- 
selection sector) in the sense that they share the same 
braiding properties: 

l'~l"4^1'-l" = Kl, l,l',l"eZ^. 

Therefore different quasiparticle types correspond to in- 
equivalent integer vectors 1 G in a iV-dimensional 
lattice, where the Bravais lattice primitive vectors are 
nothing but the column vectors of matrix K. As a re- 
sult |detK|, the volume of the primitive cell in 1-space, 
counts the number of different quasiparticle types (or su- 
perselection sectors) in a topologically ordered system 
described by effective theory (jlj. 



B. Edge excitations of an Abelian topological order 



There is a bulk-edge correspondencd^^^ for effective 
theory (iJ. When put on an open manifold Ai with 
a boundary dAi, the gauge invariance of effective La- 
grangian ^ implies the existence of edge states on 
the boundary dAi. The A^ chiral boson fields {4>i ~ 
(j)i + 27t\1 < I < N} capture the edge excitations. To 
be specific, assuming the manifold Ai covers the lower 
half-plane r2 < 0, then edge excitations localized on the 
boundary dA4 = {(fi,f2)|f'2 = 0} has the following ef- 
fective Lagrangian 



19 J 



(5) 



r2 > 0, the corresponding edge theory becomes 
^iE^-]-J2 i^i.jdohdicl^j + Vjjd,cj)id,^j). (6) 



I, J 



where IE means left edge here. V is a positive-definite 
real symmetric N x N matrix, determined by microscopic 
details of the system. The edge effective theories ([5|-([6]) 
imply the following Kac-Moody algebraP^ of chiral boson 
fields: 



r{x),dy<j>j{y)]=±2TTKjyS{x-y). 



(7) 



where +(— ) sign corresponds to the right(left) edge. The 
signature (ri+ , n_ ) of matrix K now has a clear physical 
meaning from ([5|-([6]): each positive (negative) eigenvalue 
of K corresponds to a right-mover (left-mover) on the 
right edge (|5| and a left-mover (right-mover) on the left 
edge ([6|. 

Similar to the quasiparticle excitations in the bulk la- 
beled by their gauge charge 1, associated quasiparticles 
on the edge V\ = exp{iJ2i h4'i) £^re also labeled by an 
integer vector 1 = {li, ■ ■ ■ ,/jv)"^- This identification be- 
tween bulk quasiparticle 1 and edge excitations Vi indi- 
cates that each (local) microscopic degree of freedom Q 
in the bulk also has a correspondent local excitation on 
the edge: Vj = Vki- For a, N x N matrix K, all these lo- 
cal excitations on the edge are composed of the following 
N independent local excitations (microscopic degrees of 
freedom on the edge): 

e' 51;./ K/,. 70.7(2:, t)^ 1 < / < A^. 

In the context of fractional quantum Hall states, these 
local operators on the edge are called"^^ '^^ "electron oper- 
ators". 

Here let's go over the simplest case with no symme- 
try, when symmetry group Gs = {e} and e denotes the 
identity element of a group. In this case all the (local) mi- 
croscopic boson degrees of freedom can condense in the 
bulk, and accordingly on the edge the following Higgs 
terms can be added to Lagrangian (|5|-([6]) 

Cmgas^j:iCi{e'^'Mi+ h.c.) 

= Ef=i cos {pi J2j ^LJ(t>j{x, t) + xi) ■ (8) 

where C/ and xi ^re all real parameters. Notice that 
constant factor 

PI = {3- {-l)^' ')/2, yi<i<N. 

guarantees the self statistics Q of local quasiparticle 

Mi{x,t) = e'P'^-'^''-"f'-''^'''^\ 1<I<N. (9) 



where rE stands for the right edge. On the other hand, 
if the manifold A4 instead covers the upper half-plane 



is bosonic, since if M/ is fermionic the Higgs term (|8| will 
violate locality. The Abelian topological order (featured 
by GSD on genus-g Riemann surfaces and anyon statis- 
tics) will not be affected by these Higgs terms^^ '^^'^ . 
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since all anyon excitations are local with respect to the 
microscopic boson degrees of freedom. As a result the 
condensation of local bosonic degrees of freedom {M/} 
will not trigger a phase transition, when there is no sym- 
metry in the Abelian topological order. Hence in a gen- 
eral ground these Higgs terms (8 1 should be include in 
the low-energy effective theory ( 1 1,(|5])-(|6| of an Abelian 
topological order, in the absence of any symmetry. 

Since Higgs terms ([8| are generally present in the edge 
effective theory (when there is no symmetry), they will 
introduce backscattering processes on the edge. A natu- 
ral question is the stability of gapless edge excitationi^. 
When n+ ^ n_ for the signature n_) of matrix K, 
there is a net chirality for edge states ([5|-([6| and they 
cannot be fully gapped out by the Higgs terms ([s]). A 
physical consequence is a nonzero thermal Hall conduc- 
tance in the systerrP^. If = n_ on the other hand, 
there is no net chirality on the edge. But this doesn't 
mean the edge states can be gapped out by Higgs term 

( 8 1 : the simplest counterexample is K , whose 

edge cannot be gapped out even in the absence of any 
symmetry. When the system preserves symmetry Gs, 
the structure of edge states is richer. Typically some 
Higgs terms in (|8| will be forbidden by symmet ry, and 
there will be symmetry-protected edge excitation^^sMHi] 
in the Abelian topological order. In other words certain 
branches of edge excitations will either remain gapless 
when symmetry Gs is preserved, or become gapped out 
when symmetry Gs is spontaneously broken on the edge. 
For a general discussion on the stability of edge modes 
in an Abelian topological order we refer the readers to 
section HI of Ref. [551 For the SET phases studied in this 
work, their edge stabilities are briefly discussed in section 
EHFl 



C. Different Chern-Simons theories can describe 
the same topological order 

For symmetric unimodular K matrix with det K = ±1, 
the ground state of system ([ij is unique on any closed 
manifold. Consistent with the nondegenerate ground 
state on torus, any quasiparticle 1 is either bosonic or 
fermionic with trivial mutual statistics with each o ther. 
Hence there is no topological order in the systenJSHH 
when detK = ±1. However the corresponding gapped 
phase can still have gapless chiral edge modes on its 
boundary, which are stable against any perturbations. 
Well-known examples are the integer quantum Hall ef- 
fects where K is an x identity matrix. On the other 
hand, if 

Trivial phase : dctK = (-l)^/^ = dimK = evei(10) 

the edge excitations will be non-chiral (the same number 
of right- and left-movers) and are generally gapped in 
the absence of symmetrj^. In these cases we call the 
corresponding phase a trivial phase in 2+1-D, since it's 



featureless both in the bulk and on the edge and can 
be continuously connected to a trivial product state 
without any phase transitioiJ^. 

One key point we want to emphasize is that the Chern- 
Simons theory description for a certain topologically or- 
dered phase is not unique. In other words, two different 
K matrices for effective theory ([T]) can correspond to the 
same topological phase, with the same set of quasiparti- 
cle (anyon) excitations. The two features described be- 
low are crucial for the classification of symmetry enriched 
topological orders. 

First of all, the following GL{N,Z) transformation on 
the K matrix yields an equivalent description for the 
same phase 



K ~ K = X^KX, X e GL{N, Z). 



(11) 



where GL{N,Z) represents the group of N x N uni- 
modular matrices. This GL{N,'L) transformation X 
merely relabels the quasiparticle (anyon) excitations so 
that 1 -> 1 = X~^l. It's straightforward to see that all 
the topological properties, such as quasiparticle statistics 
and GSD are invariant under such a GL{N, Z) transfor- 
mation. A brief introduction to GL{N, Z) group is given 
in Appendix \K\ 



Secondly, notice that a trivial phase satisfying ( 10 1 can 
always be added to a topologically ordered phase with- 
out changing any topological properties (such as quasi- 
particle statistics, GSD and chiral central charge of edge 
excitationd^Sl). One just needs to enlarge the Hilbert 
space to include some new microscopic degrees of free- 
dom, which form a trivial phase. Mathematically addi- 
tion of a topologically ordered phase with matrix K and 
a trivial phase with matrix Kj satisfying ( 10 1 is carried 
out by the matrix direct sunP^ 

K~K'=KeKt, (12) 
detKt (-1)^*/^^ Nt = dimKf = even. 

Therefore two K matrices of different dimensions can de- 
scribe the same topologically ordered phase. Typically 
in a bosonic system (where the microscopic degrees of 
freedom ar e all bosons) the generic trivial phase is rep- 
resented bjiiisni 



K, = 



Meanwhile in a fermionic system, both (|13|) and 



1 

-1 



1 
-1 



1 

-1 



(13) 



(14) 



together represent a generic trivial phase. 



D. Implementing symmetries in Abelian 
topological orders 

Our discussions in the previous section didn't as- 
sume any symmetrj!^in the topologically ordered phase. 



7 



Without any symmetry, an Abelian topological order is 
fully characterized by its K matrix. In the presence of 
symmetry, however, K matrix alone is not enough to 
describe a symmetry enriched topological (SET) phase: 
e.g. distinct SET phases that are separated from each 
other by phase transitions can share the same K ma- 
trix. The missing information is how the bulk quasi- 
particles (with currents jj) in effective theory (jlj trans- 
form under the symmetry. The corresponding informa- 
tion in the edge states ([5])-([6]) is how the chiral boson 
fields {4>i, 1 < / < N} transform under symmetry. 

We will restrict to unitary and anti-unitary onsite sym- 
metries in this work. By onsite symmetries we mean the 
local Hilbert space is mapped to itselP^ under the sym- 
metry transformation, so that the symmetries act in a 
"on-site" fashion. In this case studying the symmetry 
transformations of bulk quasiparticles (with currents jj) 
is equivalent td^ studying the symmetry transformations 
of edge chiral bosons {<j)j, I < I < N}. Henceforth we'll 
focus on the chiral boson variables on the edge to study 
their transformation rules under symmetry operations, in 
the presence of a symmetry group Gs. 

Most generally, under the operation of symmetry group 
element g ^ Gs, the chiral boson fields {0/} transform 
in the following wajB^l 



0,(x,t)^E,7'7^W?j(/)j(a:,t) 



(15) 



r;9K = (Wf) KW9, W9 e G'i(A^,Z). 



where rj^ = +1(— 1) for a unitary (anti-unitary) on-site 
symmetry. This is simply because under an anti-unitary 
symmetry operation (such as time reversal t — >■ —t) the 
Chern-Simons term e^'^^a^^d^ai changes sign, and in or- 
der to keep the Lagrangian ([ij in the bulk or ([5|-([6]) 
on the edge invariant, K must change sign under the 
GL{N,Z) rotation W^. 

Notice that the above symmetry transformations 
{W^ ,S(p^\g g Gs} must be compatible with group 
structure of symmetry group Gs . This provides a strong 
constraint on the allowed choices of GL{N,Z) rotations 
{W9} and U{1) phase shifts {(5(^f ~ (5(/)f + 2tt}. To be 
precise, the consistent conditions for symmetry transfor- 
mations {W^ ,6(p^\g e Gs} on an Abelian topological 
order characterized by matrix K is summarized in the 
following statement: 

The (nonlocal) quasiparticle excitations {Qi(x,t) = 
gi07(2;,t)| transform projectively under symmetry group 
Gs, while the (local) microscopic boson degrees of 
freedom {Mi{x,t) = e'P'^ "I' ^*'>} in (jg) must 
form a linear representation of symmetry group Gs. 
Here constant factor p/ = 1 if K/ / = even, or = 2 if 
Kij = odd. 

In the following we'll discuss why (local) microscopic 
boson degrees of freedom must form a linear representa- 
tion of symmetry group Gg . Imagine an Abelian topolgi- 
cal ordered phase preserves symmetry Gs . For simplicity 



let's consider Gs = Z2 = {g, e} for an illustration. We 
denote the generator of the Z2 group as g. It satisfies 
the following Z2 multiplication rule: 



(16) 



And under this Z2 symmetry operation g the edge chi- 
ral bosons transform as (15). Consider we weakly break 



the Z2 symmetry without closing the bulk energy gap (no 
phase transition) . Now Z2 operation g is not a symmetry 
anymore and there is no symmetry in the system. There- 
fore all the local bosonic degrees of freedom {Mj{x, t) = 
^ipiT.jiii,.i4>ji^,t)\i < I < N} can condensed and Higgs 
term (l8| should be allowed. At the same time, notice 
that g = e is still a "symmetry" of the system. When 
symmetry operation g act twice, its transformations on 
chiral bosons (/)(a;, t) = (0i(a;, t), ■ ■ ■ , 4'n{x, i)) become 

${x, t) ri9W9${x, t) + 5$B 

(W9)V(a;,i)+(liVxAr+?7SW9)<5^9. (17) 



where In-kn denotes an N x N identity matrix. And 
we must require all Higgs terms ([s]) with arbitrary pa- 
rameters {Ci,xi} are allowed by "symmetry" g^ = e in 
( 17 1. In other words all the Higgs terms in ([s]) should re- 
main invariant when Z2 operation g act twice as in (17l! 



This means the argument of any cosine (Higgs) terms in 
([s]) must be invariant up to a 27r phase, leading to the 
following relation: 



PK(Wsr=PK, PK(1 



NxN 



r]9W9)S(l)9 = 2™. 



where we defined NxN diagonal matrix P/,j = piSi^j 
and n = (ni,- • • ,njv)"^ S is an integer vector. The 
above relation can be rewritten as 

(Waf = 1^^^, 7j9K = (W9)'^KW9, (18) 
(Ijvxw + ?7»W9)J^9 = 27r(PK)-in, n e Z^. 

These are the group compatibility conditions on the sym- 
metry transformation (15 1 for a Z2 symmetry group 
Gs — {g,e = g^}. These conditions will be applied in 
the examples later. 

In a generic case, symmetry group Gs (and its multi- 
plication table) is fully determined by a set of algebraic 
relations 



— f/l 4/2 gN„ 



where {^i, • • • ,gNg} a set of generators in group Gg. 
Each algebraic relation Ami,--- ,mN gives rise to a con- 
sistent condition with an integer vector ... , just 
like (18) in the Gs = Z2 case. When all these group 



compatibility conditions are satisfied, any local bosonic 
degrees of freedom 

Bi(a;,i) =exp(iFPK(^(a;,t)), 1 e Z^. 
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is invariant under symmetry operation Ami,--- ,mN ■ By 
definition they form a linear representation of the symme- 
try group Gg. On the other hand, a generic quasiparticle 
excitation 



Vi{x,t) = exp (^H'^(f>{x,t) 



1 G Z 



TV 



could still transform nontrivially under consecutive sym- 
metry operation Ami,--- ,mNg (which equals identity e in 
symmetry group Gg). Therefore these (fermionic or any- 
onic) excitations transform projectivelj^^^HSI under sym- 
metry group Gs- 



E. Criteria for different symmetry emiched 
topological orders 

In the previous section we discussed the consistent con- 
ditions on the symmetry transformations on the quasi- 
particle excitations in an Abelian topological order. 
In terms of chiral boson fields (/)(x, t) which captures 
the quasiparticle contents in an Abelian topological or- 
der, under symmetry transformations ( |15[ ) (labeled by 
{W^ ,Sct)^\g G Gs} for symmetry group Gs), the (local) 
bosonic degrees of freedom (|9| transform linearly while 
(nonlocal) anyonic degrees of freedom {e"^^} can trans- 
form projectively. Together with matrix K which con- 
tains all the topological properties, the following set of 
data 



[K,{7^3,W^,6cf,3\geGs}] 



(19) 



fully characterizes a symmetry enriched topological 
(SET) phase in the presence of symmetry group Gs- 

A natural question is: is such a data a unique finger- 
print for a SET phase? Can two different sets of data 
describe the same SET phase? Not surprisingly the an- 
swer is yes. A trivial example is discussed earlier when 
symmetry group is trivial Gs = {e} and two different 
K matrices corresponds to the same Abelian topological 
order. So how can we tell whether two sets of data ( 19 ) 



describe the same SET phase or not? In the following 
we'll propose a few criteria, which thoroughly address 
this question. 

The first criterion comes from the physical picture that 
there is no "smooth" boundary condition under which we 
can sew two different SET phases with the same topo- 
logical order and symmetry group Gs- This is rooted in 
the fact that two different SET phases cannot be continu- 
ously (no phase transitions in between) connected to each 
other without breaking the symmetry. The above physi- 
cal picture can be made more precise mathematically in 
the following way. First we require ~ describe 
the same Abelian topological order in the absence of sym- 
metry, i.e. they have the same topological properties such 
as GSD (jdetK^I — |detK-^|) and quasiparticle statis- 
tics. This is because two SET phases are certainly dif- 
ferent if they correspond to different topological orders 
when symmetry is broken. 



Consider a left edge ^ of SET phase 
[K^,{r]3,Wl,S(fl\g G GJ] and the right edge 
of SET phase [K^, {ryS, , 5<^» |g G Gs}] are sewed 
together by introducing tunneling terms between the 
two edges (see FIG. [ij. We denote the chiral boson 
fields as {(f>j} on the left edge and {(f>j^} on the right 
edge. Notice that only microscopic degrees of freedom 
^ whose mutual statistics ([3| with any quasiparticle 
are multiples of 2tt, can appear in the tunneling term 
between the right and left edges'^ as shown in FIG. [l] 
Therefore a general tunneling term has the following 
Lagrangian density 

-Htunnel = Ea To. COS ((1^)^K^0^ - (1^)^K«0« + 

It 1« ^ 0, (l^)^K^l^ - (l^)^K«lf = 0, V a. (20) 

where Ta,ipa are real parameters. According to Kac- 
Moody algebra ([t]) for the chiral bosons, the condition 
on the 2nd line means the variables in each cosine term 



of ( 20 1 commute with itself and can be localized at a clas- 
sical value {{l^'^K^S^ 



Of course every 



tunneling term in ( 20 1 must be allowed by symmetry, 
i.e. t hey remain invariant under symmetry transforma- 
tion ([l5| . The edge states is fully gappecP^ if each chiral 



boson field 



is either pinned at a classical value or 



doesn't commute with at least one variable of the cosine 



terms in ( 20 ) . Notice that each cosine term in ( 20 1 must 



contain local operators from both edges. 

Let's take a look at the simplest case, when the two 
SET states share exactly the same set of data (19 1. In 



this case the tunneling term ( 20 1 essentially sews the left 
and right edge of the same SET phase. The smooth 
sewing between left and right edge basically tunnels the 
same local microscopic degrees of freedom Vj^j of one 
edge with its counterpart Vj^j on the other edge. The 
following tunneling term 



N 

tunnel ^ ^ 
1=1 



r,cos(^K,„,((/.^-<^«) + ^,). (21) 



is allowed by symmetry Gs and will gap out the edge 
states. Notice that all cosine terms commute with each 
other, so they can be minimized simultaneously. One 
important feature of the above tunneling terms is that 
there are | det_R'| inequivalent classical minima**^ for the 
{(pj — (j)j'} variables of the cosine terms. In other words, 
the chiral bosons will be pinned at one of the | det K\ clas- 
sical values by the above tunneling terms. These | det K\ 
have a one-to-one correspondence to the |detisr| degen- 
erate ground states on torus of the Abelian topological 
order here. 

When the two sets of data, [K^ , {7f ,Wl, 6$l\g G 
Gs}] for the left edge and [K^,{ria,W%6$%]g G GJ] 
for the right edge correspond to two different SET 
phases, on the other hand, there is no way to smoothly 
sew the left and right edges together. In this case 



when tunneling term ( 20 1 allowed by symmetry Gs 
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is added, either certain chiral boson modes remain 
gapless or the number of classical minima is more than 
I detK^l = I detK^I. And our first criterion is 



Criterion I: Two sets of data [K^ 



,{7^s,Wl5<j>l\ge 



Gs}] (for the left edge) and [K^, {yyS, , G GJ] 

(for the right edge) belong to the same SET phase if and 
only if there exists a tunneling term (20) connecting the 



two edges, which gaps out all chiral boson fields and has 
I det K^/^ I -/oZd degenerate classical minima. 

This criterion applies universally to both unitary and 
anti-unitary symmetries (such as time reversal symme- 
try). When detK^/^ = ±1 it automatically reduces to 
the criterion for different symmetry protected topolog- 
ical (SPT) phases in the Chern-Simons approacfP^. A 
direct consequence of Criterion I are the following two 
corollaries Corollary I: 

// the two sets of data share the same matrix 
= K^, and all their local microscopic degrees of 
freedom ^ transform in the same way under symmetry 
Gs, then they belong to the same SET phase since their 
edges can be sewed together smoothly by term [21). 



Next, notice that a GL{N,Z) transformation (11) can 
always be performed on a K matrix without changing 
the topological order. It simply relabels different quasi- 
particles. Besides, U{1) gauge transformations can al- 
ways be performed on gauge fields and chiral bosons 
{4>i}- The most general gauge transformations that re- 
label quasiparticles have the following form 

t)^J2 t) + A(bi, X e GLiN, Z)(22) 



where A(f>j e [0, 27r) are constants. We denote such a 
gauge transformation as {X, A^}. Under such a gauge 
transformation, the set of data (Il9]) changes as 



symmetry Gs) to a gauge fields (with gauge group 
Gs). One conjecture is that different SPT phases 
with Gg symmetry always leads to distinct Gg gauge 
theories. Naively one would expect this also holds for 
SET phases -'^^i two different G^-symmetry-enriched Gg 
topological orders (gauge theories) will lead to different 
G topological orders (gauge theories) where G/Gg = Gs, 
when the symmetry Gs is gauged. However as will be 
shown by the examples studied in this work, this is not 
true. In spite of that, the converse (inverse) statement 
is still true: 

Criterion II: After gauging the unitary symmetry 
Gs, if two SET phases (with symmetry Gs) lead to two 
different topological orders, they must belong to two 
distinct SET phases. 

Criterion II only applies to unitary symmetries. As 
will become clear in the examples, in certain cases 
(which we call "unconventional" SET phases), gauging 
an Abelian symmetry in an Abelian topological order 
will lead to non- Abelian topological orders^^. 



In the following sections we will demonstrate these 
criteria by classifying different SET phases with (anti- 
) unitary Z2 symmetries. The Abelian topological orders 
that will be studied include Z2 spin liquid^^ with K ~ 

2 Q j , double semion theorjEHMl with K ~ ( ^ ^ 

and Laughlin states'^ (K ~ to) at different filling frac- 
tions i> — 1/to. Among them Z2 spin liquid and dou- 
ble semion theory are non-chiral Abelian phases, in the 
sense that their edge excitations have no net chirality. 
And in the absence of symmetry their edge excitations 
will generically be gapped. On the other had Laughlin 
states are chiral Abelian phases with quantized thermal 
Hall conductance. 



ygeGs, W9 



K ^"^^^ X^KX 

{X,A<?} 



(23) 



X-1W9X, 



and r]^ remains invariant. Here comes the second 
corollary 



Corollary II: any two sets of data (19) that can be 
related to each other by a gauge transformation (23) 
correspond to the same SET phase. 

Last but not least, an important lesson from studies 
of SPT phases is that there is a dualitjESESi between 
SPT phases and gauge theories (or intrinsic topological 
orders). This duality is established by gauging the 
(unitary) symmetry Gs in the SPT phase, i.e. coupling 
the physical degrees of freedom (which transform under 



III. EXAMPLES 

In this section we'll apply the Chern-Simons approach 
discussed in previous sections to various Abelian topolog- 
ical orders. We start by classifying Z2 spin liquid with 
time reversal symmetry Gs — Z2 and with a unitary Z2 
symmetry Gs = Z2. Usually by Z2 spin liquids people as- 
sume that there are certain fractionalized quasiparticles 
carrying spin quantum numbers, coined "spinous" and 
other fractionalized quasiparticles carrying no spin quan- 
tum numbers, coined "visons". The mutual (braiding) 
statistics of a spinon and a vison is semionic {6s,v = "■), 
while the self statistics of a vison is bosonic. A Z2 spin 
liquid has 4-fold GSD on a torus. All these topological 
properties are captured by the Chern-Simons theory ([T]) 
with 



K 



2 
2 



(24) 
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In the context of this work, we don't assume spin ro- 
tational symmetry and vision/spinon generically cannot 
be distinguished from their quantum numbers. Despite 
this fact we still use the name "Z2 spin liquid" to label 
this Abelian topological order. The 4 degenerate ground 
states on a torus correspond to the 4 superselection sec- 
tors, which are associated with the 4 inequivalent quasi- 
particles: 



(25) 



where both e and m have bosonic self statistics and 
they denote electric and magnetic charge in a Z2 gauge 
theorjEH respectively. / is the bound state of an elec- 
tric charge and a magnetic charge, with fermionic self 
statistic. corresponds to any local microscopic degrees 
of freedom Q, belonging to the vacuum sector. In the 
folklore of Z2 spin liquid, a vison is e (or m), and accord- 
ingly a bosonic spinon is m (or e). 



Classifying Z2 spin liquids with time reversal 
symmetry 

As a warmup we consider Abelian topological order 

g^} with alge- 



(24 1 with symmetry group = {g,e 
bra (16). Notice that the generator of Zj group, g is an 
anti-unitary operation with 7/^ = — 1 in ( 15 1. In this case 
we rely on Criterion I and its corollaries to differentiate 
different Z^^-SET phases. 



The associated group compatibility condition ( 18 1 for 
Gs = Zj in Abelian topological order ( 24 ) is 



(WS)' = l2x2, - 

(12X2-W9)<5<?9 



2 

2 



: TT 



1 

1 



n, neZ2. (26) 



The e GL{2,Z) solution to the above conditions is 
'1 
.0 -1, 



± 



However notice that the following 



GI/(2,Z) gauge transformations (23 1 keep the K matrix 
(I24I invariant: 



±1 



2x2, 



± 



1 

1 



Therefore Wf 



} \) - w ^ (-0' ; 

are equivale nt, related by gauge transformation X 
^0 1' 
1 



± ( " ^ I in ( 23 1 . And we can fix the gauge by choosing 



W9 = 



1 
-1 



K 



2 
2 



with symmetry Gs ~ Z2 ~ {g,e — g^} 



Data set in 



19 



2 
2 



-1,W9,<503}] 



Label 


#1 


#2 




(i " ) 


il 




(0,0)^ 




Proj. Sym. (m = (0, 1)^) 


No 


Yes 


Gapless edges 


No 


No 



Table I: Classification of difi^erent Z2 spin liquids with (anti- 
unitary) time reversal symmetry. There are two diflerent SET 
phases: in SET phase #1 all quasiparticles in 1 25 1 transform 



linearly under Z2 symmetry, while in SET phase #2 quasi- 
particle m transforms projectively under Z2 symmetry. The 
data set in the 2nd line completely characterize these SET 
phases. "Proj. Sym." is short for "projective realization of 
symmetry" in the table. 



for time reversal operation g. Then solving the 2nd line 
of conditions ( 26 1 we obtain 



n2 = 0, 



mod 2tt. 



We can alway s c hoose a gauge transformation {X = 
l2x2,A(/)} in (23) so that 6411 = 0. Meanwhile since 
1 = (2, 0)-^ and 1 — (0, 2)-^ are the local excitations in 
the system, according to Corollary II, ni =even all cor- 
responds to the same SET phase. Meanwhile ni =odd 
leads to another SET phase, which is distinctive from 
the Til =even SET phase. This is because the magnetic 
charge m ~ (0, 1)"^ transforms projectively in ni =odd 
phase but transforms linearly in the 712 =even phase un- 
der time reversal symmetry. It's straightforward to check 
that there is no way to smoothly sew the two edges of 
ni =even and ni =odd SET phases by a time-reversal- 
invariant tunneling term (20), which has 4- fold degener- 
ate classical minima. Therefore according to Criterion I 
they belong to two different SET phases. These are the 
only two different classes of 
spin liquids, as summarized in TABLE [l] 

Since the previous calculations are based on 2 K ma- 
trix (24), it is natural to ask: what if we enlarge the di- 



Z2 -symmetry-enriched Z2 



mension of K matrix ( 12 1 by introducing the trivial part 
with ( 13 1? In this new representation of the same Abelian 
topological order, will we get more SET phases or not? 
Notice that the trivial part (ll Sl is nothing but the K 

lai^I^r 



matrix for a bosonic SPT phase 



in 2+1-D. For anti- 



unitaryZj symmetry, there is no nontrivial bosonic SPT 
phas(P0|26J in 2+1-D. This means the edge chiral bosons 
for the trivial parts can always be gapped out by intro- 
ducing symmetry-allowed backscattering cosine (Higgs) 
terms, whose classical minima is pinned at a unique clas- 
sical value since |detKt| = 1. According to Criterion 
I, in the presence of Zj symmetry, when the dimension 
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under time reversal operation g. As a result the edges 
can be completely gapped by introducing Higgs terms 

'Hhiggs = CCOS {2(f)i{x , t)) . 

which pins chiral boson field (j)i{x,t) to a classical value 
{4>i{x,t)) = or TT, without breaking the time reversal 
symmetry. Therefore in general there are no gapless edge 
states for the two SET phases with Gs ~ ■ 

Potentially, one could conceive of a phase where both 
electric and magnetic charges transform projectively un- 
der time reversal symmetry. However, such a phase is 
only possible as the surface state of a 3D SPT phase 
with time reversal symmetry^-. The K-matrix classifi- 
cation correctly reproduces the fact that this phase is 
forbidden. 



B. Classifying Z2 spin liquids with onsite Z-i 
symmetry 

As discussed earlier, the reason why 2x2 matrix K ~ 



.2 0, 



is enough to describe all Z2 -symmetric Zi spin 



liquids is that there is no nontrivial Zj-SPT phases of 
bosons in 2+1-D. In other words the possible trivial part 
( 13 1 that can be added to K in ( 12 1 doesn't contribute 



to new structure to SET phases. However, for a unitary 
Gs = Z2 symmetry, as will b ecome cl ear later, there is a 
nontrivial bosonic SPT phas^^^l^^l^ whose edge cannot 
be gapped without breaking the Z2 symmetry. This Z^,- 
SPT phase can be understood in Chern-Simons approach 

^0 1^ 

,1 0, 



with K 



Therefore we need to consider 4x4 



matrix 



K = 




/O 2 0\ 
2 
1 

yo 1 oy 



to represent a generic (unitary) Z2-symmetry-enriched 
Z2 spin liquid. The group compatibility condition (18) 
for the unitary Z^ symmetry (jf = 1) transformation 



of K is enlarged by adding the trivial parts, it will not ( 15 1 becomes 
introduce any new SET phases. 

At the end we discuss the stability of edge excitations 
in the two SET phases. Notice that the chiral bosons 
{01,2} transform as 



(W9)' = l4x4, 



(27) 



20/ ll 0/ ^ ' \2 0/ \1 0/ 




The gauge inequivalent solutions to is the followingZSl 

W9 = l4x4, ^ © 12X2. 

We will discuss these two cases separately in the follow- 
ing. 

1. "Conventional" Zi-symmetry-enricheA Z2 spin liquids 



First we discuss the solution W^* 



I4 



In this 



case the anyon quasiparticles ( 25 ) in the Z2 spin liq- 
uids merely obtains a C/(l) phase factor under the Z2 
symmetry operation g, and we call them "conventional" 
SET phases. Due to the gauge transformations X = 




2x2, 1-2x2 




which leave the K matrix 



invariant, we know that the integer vector n in (27 1 has 
the following equivalency relation 



/"A 

"-2 



ni 



Moreover Corollary II tells us 



An + 






( ^ 




( ^ 




( ^ 


"2 






"2 +2 




"2 




"2 


"3 






"-3 




713+2 




"3 


V "-4 


) 




\ Hi J 




I "4 ) 




\n4 + 2j 



and 



\ 








"2 






1 





/ 




VI/ 



ni 

n2 




As a result we obtain 6 inequivalent SET phases with 
W^' = l4x4 under Z2 symmetry, with their symmetry 
transformations 5(t)^ summarized in TABLE Lll We re- 
quire 5^ ^ so that the local excitations ^ form a faith- 
ful representatiorPfil of symmetry group Gg. In the fol- 
lowing we briefly discuss the consequence of gauging the 
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with unitary symmetry G s = Z2 ~ {g , e — g^} 
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,{r;« = +1,^ = 14x4,50"}] 


Label 


#1 


#2 


#3 


#4 


#5 


#6 




(0,0,^,0)^ 


(0,0,7r,7r)^ 


(V2, 0,^,0)^ 


(7r/2,0,7r,7r)^ 


(V2,^/2,^,0)^ 


(7r/2,7r/2,vr,7r)^ 


Proj. Sym. (e ~ (1,0,0,0) ) 


iNo 


No 


Yes 


Yes 


Yes 


Yes 


Proj. Sym. (m~ (0,1,0,0) ) 


iNo 


No 


iNo 


INo 


Yes 


Yes 


Proj. bym. (/ ~ (1,1,0,0) ) 


iNo 


No 


Yes 


Yes 


INo 


INo 


Symmetry protected edge 


iNo 


Yes 


INo 


Yes 


Yes 


Yes 


After gauging symmetry: 


/O 2 0\ 
2 
2 

\0 2 0/ 


/O 2 \ 
2 
2 

\0 -2/ 




(::) 


(::) 


(: 


(: ) 


hqg mod 1 





1/4 





1/4 


1/8 


3/8 


Oqg,e mod 2tt 








7r/2 


7r/2 


7r/2 


7r/2 


dqg.m mod 2-K 














7r/2 


7r/2 


Oqgj mod 27r 








7r/2 


7r/2 




TT 


Comparison to Ref. [17] 


(000) 


(100) 


(010)&mi = 


(110) 


mi = 2 


? 



Table II: Classification of "conventional" Z2 spin liquids enriciied by onsite (unitary) Gs = Z2 symmetry. There are 6 different 
"conventional" SET phases, where under Z2 symmetry all quasiparticles (e,m, /) merely obtain a (7(1) phase factor. The data 
set in the 2nd line completely characterizes these SET phases. Kg denotes the topological order, which is obtained by gauging 
the unitary Ga = Z2 symmetry in the Z2 spin liquid. Some of these SET phases have Z2 symmetry protected edge states, which 
will be gapless unless Z2 symmetry is spontaneously broken. On gauging the Z2 symmetry (blue entries), a new quasiparticle 
excitation {Z2 vortex) qg is obtained, as described in Appendix[B] Its statistics (B3l-(B4| are also summarized in the table: its 
self statistics Oq^ = 2-nhqg has a one-to-one correspondence with its topological spin exp(27ri/iqg). 



unitary Z2 symmetry. A detailed prescription of gaug- 
ing a unitary symmetry in the Chern-Simons approach 
is given in Appendix ^ where we've shown that gauging 
Z2 symmetry {W^ = Uy.i,5$B = ^(^1/2,^2/2,1,^4)^} 
yields an Abelian topological order described by matrix 
Kg in (B2|. Take #6 for an example, after gauging 
symmetry we have 



/O 2 \ 

2 -1 
2 -1 
\2 -1 -1 -2/ 

where the first equivalency 



formation ( 23 1 with 



/4 0\ 

0-400 
1 

Vo 1 0/ 

is realized by gauge trans- 




X = 



/ 2 -1 

1 -1 

-1 3 1 

V 2 -2 





1 



detX = 1. 



From TABLE |lT] one can see that two different SET 
phases can lead to the same (intrinsic) topological or- 
der by gauging their Z2 symmetry, such as #3 and #4, 
or #5 and #6. 



The stability of edge excitations is also summarized in 
TABLE |TT] For SET phases #2 and #4 the gapless edge 

excitations come from the trivial ^\ part (lower 2x2 

part) of K matrix, which corresponds to the symmetry 
protected edge modes of bosonic Z2-SPT phases. How- 

^0 2\ 



ever for #5 and #6, the topologically ordered 



2 0, 



part (upper 2 x 2) of K matrix also contribute to c = 1 
gapless edge states. In other words in a Z2 spin liquids, 
if both e and m transform projectively under Z2 symme- 
try, the edge excitations is protected to be gapless un- 
less symmetry is broken. The edge chiral bosons {0i.2} 

^0 2^ 



for K 



,2 0, 



can be refermionized as right-moving 



branch ipR ~ exp [i(0i + (^2)] and left-moving branch 
tAl ~ exp [i(0i — 02 )] ■ The edge effective theory ([H]) can 
be rewritten as 

CrE = - v+dx)ijR ~ iiplidt - v^dx^JL- (28) 



where v± = (Vi^2 ± ^Vf_2 + Vi^iV2,2)/2 are the ve- 
locities of edge modes. It's easy to see that under Z2 
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2 
2 



nth unitary symmetry Gs = Z2 ~ {g, e — g^} 



2 
2 



1 

1 



+1, w 



1 

1 



Label 



#1 



#2 



#3 



#4 



(0,0,^,0)^ 



(0,0,7r,7r) 



(7r/2,V2,^,0)^ 



(7r/2,7r/2,7r,7r)^ 



Proj. Sym. (e : 














1 

















No 



No 



Yes 



Yes 



Proj. Sym. (/ ~ (1, 1, 0, 0)^ 



No 



No 



No 



No 



Symmetry protected edge states 



Yes 



Yes 



Yes 



Yes 



Central charge (c mod 1) of gap less edge states 



1/2 



1/2 



1/2 



1/2 



ha„ mod 1 



16 ' 16 



16 ' 16 



16 ' 16 



Relation to Kitaev's 16-fold waj** 



{v = l)®{v= 15) 



{v = ^)®{v= 11) 



(jy = 7) (g) (i/ = 9) 



(i/ = 3) (g) (1/ = 13) 



Table III: Classification of "unconventional" Z2 spin liquids enriched by onsite (unitary) Gs = Z2 symmetry. There are 4 
different "conventional" SET phases, where under Z2 symmetry quasiparticles e and m will exchange. The data set in the 2nd 
line completely characterizes these SET phases. Kg denotes the topological order, which is obtained by gauging the unitary 
Gs = Z2 symmetry in the Z2 spin liquid. All these SET phases have Z2 symmetry protected edge states, which will be gapless 
unless Z2 symmetry is spontaneously broken. However, the central charge c of (symmetry protected) gapless edge states are 
different, qg is the new quasiparticle excitation (the Z2 vortex) obtained by gauging the "unconventional" Z2 symmetry in 
the system, as described in Appendix [D] Gauging this "unconventional" Z2 symmetry leads to new non-Abelian quasiparticles 
(blue entries), which has quantum dimension dq^ = y/2 and topological spin exp(27ri/!,qg ). The quasiparticle contents of the 
non-Abelian topological orders obtained by gauging Z2 symmetry is summarized in TABLE |VI| The "gauged" non-Abelian 
topological orders for all these 4 SET phases have 9-fold GSD on a torus, corresponding to 9 different superselection sectors. 



symmetry g the chiral fermions transform as -0^ — >■ —^b, 
and ijjL — ^ V'i for SET phase #5, and hence backscatter- 
ing term 4'L'4'R + h.c. and V'lV'i? + ^-c are forbidden by 
Z2 symmetry. As a result there are gapless edge states 
in ^5 SET phase, protected by Z2 symmetry. 

From TABLE |lT] it seems that in SET phases #5 
and #6, anyonic quasiparticles (e, m, /) transform in the 
same fashion. They both have symmetry protected edge 
states, and gauge into the same intrinsic Z2 topological 



order K„ 



4 0, 



It is Criterion I that dictates they 



are different SET phases. Physically they do have dif- 
ferent edge state structures: both in #5 and ^6 a pair 
of counter-propagating chiral fermions from topological 

order part are protected by symmetry, while in 

7^6 there are also symmetry protected edge modes from 

the trivial | ^ "'^1 part. 
U 0/ 



2. "Unconventional" Z2- symmetry- enriched Z2 spin liquids 



Here we focus on solutions to (27 1 with 
^0 1^ 



1 0, 



l2x2- Notice that one can always choose a 



(50f , and hence 



proper gauge A0 in (23 1 so that 
ni = n2 in (27). There are 4 inequivalent Scj)^ solutions of 



this type to (27), as summarized in TABLE III We can 



see that the electric charge e and magnetic charge m form 
a two-dimensional representation of Z2 symmetry group. 
They transform projectively in SET phases #3 and #4, 
or transform linearly in SET phases #1 and #2. All 
these phases host symmetry protected edge excitations 
on the boundary, but the structure of these gapless edge 
states (when symmetry is preserved) are different. 

For SET phase ^1 and ^3, the edge chiral bosons 

^ ^ part of K matrix can be fully 



"3,4 



for the trivial 




However, the chiral bosons 01 2 
^0 2^ 

by Z2 symmetry. To be precise, in the refermionized de- 



gapped by a cos < 
for topologically ordered 



part of K are protected 



scription (28) for edge states, the chiral fermions trans- 
form as 



0j 



4- 



(29) 



for (50f = 7r(ni/2, ni/2, 1, ng) with rii^g = 0, 1 in TABLE 
mil We can rewrite each chiral fermion in terms of two 
Majorana fermiorP^ ^^/l and t^/j/l 

For phases #1 and #2 with ni = 0, the following 
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backscattering term is allowed by Zi symmetry 
Ubs oc + V^l) + h.c. = 4^fl^L. 

Therefore the ^/j/^ branch of Majorana fermions are 
gapped, while the ruiji^ branch is protected by sym- 
metry. For phases #2 and #4 with ni = 1, similarly the 
following backscattering term 

'H.hs oc -^flX^L - ■^\) + /^•c. ^ ^mmh- 

is allowed by Z2 symmetry. It will gap out TiRjh 
modes and leave Majorana modes ^i^/^ gapless. As a 
consequence a c = 1/2 branch of Majorana fermions 
will remain gapless, unless Z2 symmetry is sponta- 
neously broken on the edge. Together with the c = 1 

Z2-symmetry-protected chiral boson edge^*" from 

part for phases #2 and #4, we obtain the total central 
charge c for all 4 "unconventional" SET phases as 
summarized in TABLE IIIII 



Aside from gapless edge states, another important fea- 
ture for these "unconventional" SET phases is that they 
lead to non-Abelian topological orders once the unitary 
Z2 symmetry is gauged. For these unconventional SET 
phases, a vertex algebra approach is introduced in Ap- 
pendix [D] to gauge the unitary symmetry. The quasipar- 
ticle contents of the "gauged" non-Abelian topological or- 
ders for these 4 SET phases are summarized in TABLE 
|VI[ The "gauged" topological orders are related to the 
"unconventional" Z2 gauge theories describing fermions 
with odd Chern number v coupled to a Z2 gauge field, 
as Kitaev described in his 16- fold way classification of 
2+1-D Z2 gauge theorie^^H, More specifically, these non- 
Abelian topological orders are Z^ x Z2 gauge theories, 
the direct product of v =odd Z^ gauge theory and its 
time reversal counterpart P = 16 — i^, as summarized in 
TABLE |III[ Hence these "gauged" topological orders all 
have non-chiral edge states (chiral central charge c_ =0), 
which will generally be gapped in the absence of extra 
symmetry. 

After gauging the symmetry, new quasiparticles with 
quantum dimension dq^ = \/2 emerge as deconfined exci- 
tations, called Zi vortices. When any quasiparticle q in 
the original SET phase is moved adiabatically around a 
Z2 vortex once, it becomes its image gq under Z2 sym- 
metry operation. These Zi vortices are similar to a Ma- 
jorana bound state in the vortex core of a spinless p + ip 
superconductor'^ in 2+1-D. However, they have different 
topological spin than those in ip superconductors. To 
be specific, there are 4 inequivalent Z^ vortices with topo- 
logical spin exp(±Y^27ri) and exp(±^^27ri), as shown 
in TABLE |III[ All these non-Abelian topological orders 
have 9-fold GSD on a torus, corresponding to 9 different 



superselection sectors shown in TABLE VI It's not hard 
to see that SET phases #1 and #3 lead to the same 
non-Abelian theory by gauging the Z^ symmetry, and 
so do SET phases #2 and #4. However, SET phases 



#1 and #2 do lead to different non-Abelian topologi- 
cal orders after Z^ symmetry is gauged. In particular, 
their 9x9 modular S and T matrices in the basis of 
TABLE |VT] are shown in the end of Appendix |D] Com- 
bining the edge states and "gauged" topological orders 
summarized in TABLE |III[ indeed there are 4 distinct 
"unconventional" Z2-symmetry-enriched Z2 spin liquids. 



3. Discussions on Z2-symmetry-enriched Z2 gauge theories 

First we discuss the results for "conventional" SET 
phases with onsite Z2 symmetry. For Z2 spin liquids 

^0 2^ 
2 0^ 

ferent conventional Z2-SET phases as summarized in TA- 

'2 ^ 

on-site Z2 symmetry we obtain 8 different conventional 
Z2-SET phases as summarized in TABLE |TV] Here we 
make some connection between these SET phases and 
Z2-symmetr y-enri ched Z2 gauge theories obtained in pre- 



K 



BLE 



with on-site Z2 symmetry we obtain 6 dif- 



II 



For double semion theory K 



with 



3tr y- 
vious studie^i^'i^. 

In Ref. [T5]the exact soluble lattice models for 8 differ- 
ent Gs = .Z^2-symmetry-enriched Gg = Z2 gauge theories 
are obtained. They correspond to the following group co- 
homology n^{Gs x Gg, U{1)) = n^{Z2 x Z2, U{1)) = Zl 
Among these 8 different SET phases, 4 comes from Z2 
spin liquids with on-site Gg = Z2 symmetry, and the 
other 4 from double semion theory with onsite Gs = Z2 
symmetry. They are nothing but #1, #2, #3, #4 in TA- 
BLE |TT] together with #1, #2, #7, #8 in TABLE [TVl 

In Ref. |T7] twelve different Z2-symmetry-enriched Z2 
topological orders are proposed, among which six are Z2 
spin liquids and the others are double semion theories. 
It was conjectured that different G^-symmetry-enriched 
Gg-gauge theory (with gauge group Gg) are classified 
by group cohomology W^'^^{G, U{1)) in d-spatial dimen- 
sions, where G is an extension of symmetry group Gs 
by gauge group Gg. When Gs ~ Gg — Z2 we have 
G/Gs = Gg and hence G = Z2 x Z2 or G Z4. The 
number 12 = 2^-1-4 is associated with ^^(^2 x Z2, J7(l))® 
n^iZi, [/(I)) = Z^©Z4. The 8 different SET phases from 
'H'^{Z2 X Z2, U{\)) are the same as in Ref. [16] which are 
discussed earlier. After gauging the Gs = Z2 symmetry, 
these 8 different SET phases lead to Abelian topological 
orders described by a 4 x 4 matrijP^ 



K(nirt2n3) 



(2ni 2 n2 0\ 

2 

n2 2^3 2 

V 2 0/ 



where ni,n2,ri^ = 0,1. It's not difficult to check that 
these 8 different SET phases labeled by (71171271.3) have 
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2 



with unitary symmetry Gs = Z2 — {g, e — g^} 
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e 



1 

1 

^3 



,{7]'' = +l,W'' = l4x4,<50''}] 

#4 #5 ^ 



U/2\ / \ 

7r/2 

V / V / 

Yes No 



#7 



#8 



Proj.Sym.(s) 





TT 

Vo/ 

No 





TT 

No 





TT 

Yes 



7r/2 

TT 

V w 
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^^/2\ 

7r/2 

TT 

Yes 



/7r/2\ 

7r/2 

TT 

Yes 



Proj.Sym.(s) 



No 



No 



No 



No 



Yes 



Yes 



Yes 



Yes 



Proj.Sym.(b) 



No 



No 



Yes 



Yes 



Yes 



Yes 



No 



No 



SP edge 



No 



Yes 



Yes 



Yes 



Yes 



Yes 



No 



Yes 



On gauging the symmetry 



hq^ mod 1 



(2 0^ 

0-200 
2 

Vo 20/ 





(2 Q Q\ 
0-200 
2 

Vo -2/ 

1/4 



8 
-2 

1/16 



8 
-2 



2 




5/16 -1/16 



2 

-i 

3/16 



4 
4 



/O 2 0\ 
2 2 10 
12 2 

Vo 2 0/ 

1/4 



mod 27r 











7r/2 



7r/2 







7r/2 



7r/2 



mod 27r 







7r/2 



7r/2 



7r/2 



7r/2 



mod 2-K 



■k/2 



■k/2 



Tx/2 



7r/2 



Notation in Ref. [TTl 



(001) 



(101) 



mi = 3 



mi = 1 



(Oil) 



(111) 



Table IV: Classification of double semion theory (CI I enriched by onsite (unitary) Gs ~ Z2 symmetry, see Appendix [C] for 
details. There are 8 different "conventional" SET phases, where under Z2 symmetry all quasiparticles (s, s, b) merely obtain a 
17(1) phase factor. The data set in the 2nd line completely characterizes these SET phases. Kg denotes the topological order, 
which is obtained by gauging the unitary Gs ~ Z2 symmetry in the double semion theory. Some of these SET phases have 
Z2 symmetry protected (SP) edge states, which will be gapless unless Z2 symmetry is spontaneously broken. On gauging the 
Z2 symmetry (blue entries) a new quasiparticle qg is obtained, as described in Appendix [B] Its statistics ([C5|-([C6| are also 
summarized in the table: its self statistics Oq^ = 2Tihqg has a one-to-one correspondence with its topological spin exp(27ri/i5g). 
Note, there are no 'unconventional' symmetry realizations with this topological order. 



the following correspondence with our results: 



and 



K(OOO) 



K(IOO) 



K(OIO) 



/^O 2 0\ 

2 

2 

VO 2 0/ 

/^O 2 \ 

2 

2 

Vo -2/ 

/o 4 



^ #1 in TABLE ini 



^ #2 in TABLE HH 



K(llO) 



A 0, 



■0 4' 
,4 0, 



^ #3 in TABLE nil 



^ #4 in TABLE HH 



K(OOl) 



K(lOl) 



K(Oll) 



K(lll) 



(2 o\ 

0-200 
2 
VO 2 0/ 

/2 \ 

0-200 

2 

VO -2/ 

^0 4 

/O 2 0\ 
2 2 10 
12 2 
Vo 2 0/ 



^#1 in TABLE IVl 



^ #2 in TABLE HV 



^ #7 in TABLE nv 



^ #8 in TABLE lia 



We want to emphasize that when the on-site unitary Z2 
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symmetry is gauged, different SET phases could result 
in the same topological order. Therefore the structure of 
the topological order obtained by gauging the symmetry 
doesn't fully characterize one SET phase! For example 
with Xi,2 e GL{4,Z) 



Xf K(010)Xi = X^K(011)X2 = 



Xi = 



/ 2 

1 

-2 

V-1 



1\ 



1 
0/ 



,Xo = 



V-1 



/O 4 0\ 

4 

1 

yo 1 0/ 

2 -1 l\ 

1 

-2 1 

1 0/ 



and it's straightforward to see K(IOO) ~ K(OOl) and 
K(llO) ~ K(Oll). In fact the 8 different matrices 
K(nin2'^3) describe 5 different Abelian topological or- 
ders. Since these 8 theories correspond to different 
Dijkgraaf-WitteiJ^theories, gauge theories with topolog- 
ical terms specified by 7i^{Z2 x Z2, U{1)) — Zf, this also 
implies that different Dijkgraaf-Witten theories based on 
a particular gauge group share the same topological or- 
der. Perhaps further information regarding which parti- 
cles comprise electric and magnetic charges is required to 
uniquely define those phases. 

The other 4 SET phases proposed in Ref . [T7] are asso- 
ciated to group cohomology T-L^{Z4,U{1)) = Z4. Ref. [T7| 
asserted that after gauging the Z2 symmetry they lead 
to Abelian Z4 topological orders described by 



K(mi 




TO = 0,1,2,3. 



We found that these 4 different SET phases have over- 
lap with those 8 SET phases associated to H^{Z2 x 
Z2, t/(l)) = they are nothing but 



#3 in TABLE HJ 
^ #5 in TABLE nil 
^ #5 in TABLE HV] 
^ #3 in TABLE nV 




"Conventional" SET phases #6 in TABLE [IT] {Z2 spin 
Hquids) and #4, #6 in TABLE [!v| however, doesn't fall 
into these 12 classes. It's not clear at this point whether 
these SET phases can be described by group cohomology 
theory. 
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Figure 2: (color online) A fermion mode (/) localized at the 
boundary between two subsystem A and B which from a bi- 
partition of the on a sphere, where the "unconventional" (type- 
II) Ising-symmetry-enriched Z2 spin liquid resides. Under the 
Ising symmetry operation, an electric charge e will transform 
into a magnetic charge m. Consider one electric charge is 
created in each subsystem. If we perform Ising {Z2) sym- 
metry only on subsystem A, a fermion mode will emerge on 
the boundary, as the electric e charge turns into a magnetic 
charge m va A. 



Measurable effects of "unconventional" Z2 symmetry 
realizations 



Suppose there is a Z2 spin liquid which preserves 
Z2 spin rotational symmetry (for integer spins), are 
there measurable effects for these SET phases? More 
specifically, what are the distinctive measurable fea- 
tures of the "unconventional" Z2-SET phases? In this 
section we'll try to answer these questions in two 
aspects, i.e. measurements in the bulk and on the edge. 
We'll focus on unconventional SET phases in this section. 

First of all, an important ingredient of SET phases is 
how their quasiparticles transform under symmetry Gg. 
This gives us one way to measure an SET phase: to cre- 
ate quasiparticle excitations and apply symmetry oper- 
ation on them. For example for a Z2 spin liquid on a 
closed manifold (a sphere, say), a pair of electric charges 
e ~ (1,0,0,0)"^ (or magnetic charges to ~ (0,1,0,0)^) 
can be created on top of the groundstate. For an on- 
site unitary symmetry (such as Z2 spin-fiip symmetry 
g), one can choose to perform the symmetry operation 
only on a part of the whole system. For example, FIG. 
|2] shows such a striking measurable effects on the un- 
conventional Z2-symmetry-enriched Z2 spin liquids. As- 
sume in the SET phase a pair of electric charges e are 
created, one in subsystem A and the other in subsystem 
B. The whole system AljB lives on a closed manifold, 
say a sphere on which the groundstate is unique. If we 
only perform the "unconventional" Z2 symmetry opera- 
tion g in subsystem A (but not in B) , the electric charge 
e ~ (1,0,0,0)"^ in subsystem A will become a magnetic 
charge to ~ (0, 1, 0, 0)"^. However one electric charge and 
one magnetic charge cannot be created simultaneously 
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a/2|±) = |e) ± |m) 




± 




Figure 3: (color online) The Ising symmetry eigenstates are 
linear combinations of minimal entropy states (MESs) for 
a "unconventional" Ising-symmetry-enriched spin liquid, 
since one MES je) transforms into another MES |m) under 
Ising symmetry operation. 



(e X m = / 7^ 1) on top of the groundstate: the conserva- 
tion of "topological charge" requires the existence of an 
extra fermion / in the system! Such a ferniion mode / 
lives on the boundary (dashed line in FIG. [2]) between 
subsystem A and i3, as shown by the wavy line in FIG. 
|2] This effect happens in all 4 SET phases in TABLE [ml 



Secondly, there are degenerate ground states once 
we put the SET phases on a closed manifold with 
nontrivial topology (with nonzero genus). For example, 
they have 4- fold GSD on a torus (or infinite cylinder). 
How these ground states transform under symmetry 
is a reflection of how anyon quasiparticles transform 
under symmetry. Specifically, one can always choose a 
set of basis called the minimal entropy states (MESsj^. 
As the name implies, each MES is a superposition of 
degenerate ground states which minimizes the bipartite 
entanglement entropj^^, once a certain entanglement 
cut is chosen: e.g. along the y-direction in the middle 
of the infinite cylinder as shown in FIG. [3] The MESs 
are fiux eigenstates^, which keeps maximum knowledge 
of the states after the entanglement cut. Specifically 
for a Z2 spin liquid, we can label the 4 MESs as 
|1), |e), |m), I/) on an infinite cylinder. Remarkably 
under the "unconventional" (type-II) Ising symmetry 
operation, two MESs (|e) and |m) exchanges) and their 
linear combinations |±) — (|e) ± |m))/-\/2 are the Ising 
symmetry eigenstates. Therefore the MESs necessarily 
breaks Ising symmetry in such an unconventional SET 
phase! This phenomena can be actually measured in 
numerical studiei^. 

Thirdly, the edge state structure encodes many in- 
formation of a SET phase, when it supports symmetry 
protected edge modes. For unconventional SET phases, 
there are always gapless edge excitations protected by 
symmetry, unless the symmetry is spontaneously bro- 
ken on the boundary. In the specific case of uncon- 
ventional Z2-symmetry-enriched spin liquids summa- 
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breaking 



Figure 4: (color online) Domain wall bound state on the edge 
of "unconventional" Ising-symmetry-enriched Z2 spin liquids 



(see TABLE III I. In these SET phases, under Z2 symmetry 
operation, one electric charge will transform into a magnetic 
charge and vice versa. The on-site unitary Z2 (Ising) sym- 
metry can be, e.g. a spin-flip symmetry. On the two sides 
of the Ising-symmetry domain wall, two different backscat- 
tering "mass" terms related by spin-flip Ising symmetry are 
added to gap out the edge states. These two mass terms 
break Z2 symmetry in opposite ways. A non-Abelian bound 
state with quantum dimension dq^ = \/2 is localized at each 
Ising domain wall. For a "conventional" Z2-SET phase, such 
a Ising mass domain wall will trap an Abelian bound state 
with quantum dimension 1. 



rized in TABLE [III| one important feature is that SET 
phases #1 and #2 supports gapless (non-chiral) Majo- 
rana fermion excitations on the boundary, with central 
charge c = 1/2 mod 1. However this is not universal 
for all unconventional SET phases. A more interesting 
effect comes from the bound state localized at a do- 
main wall on the edge. Take SET phase ^1 for example, 
a perturbation on the edge 



■Hi = A\ cos(2(/)i) -I- A2 cos(2(/!<2) + ^4 cos < 



(30) 



K = 



if Ai_4 7^ (or A2A 7^ 0). On one 



can fully gap out the edge excitations in (pi- (pi with 

side of the Z2 domain wall on the edge of SET phase 
^1, we break Zi symmetry with A2 = Q and A\a 7^ 0. 
On the other side of the Z2 domain wall, Z2 symmetry 
is broken in the opposite way so that = and 
A2A 7^ 0. Physically the electric charges are condensed 
on one side of the domain wall, while magnetic charges 
condense on the other side. At the domain wall a 
non-Abelian (Majorana) bound stat^^o] jg localized, 
which has quantum dimension v^, as illustrated in FIG. 
m Such a domain-wall-bound state is similar to those 
localized at the (ferromagnetism / superconductivity) 
mass domain walls of a quantum spin Hall insulator^il 
In the vertex algebra context, these domain- wall-bound- 
states correspond to quasiparticle qg (ii^ the 8th role) in 
TABLE IVI] In SET phase #1 e.g. it has topological spin 
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exp(— i7r/8). In the bulk-edge correspondence of SET 
phases, such a bound state on the edge is related to the 
Z2 vortex (Jo in the bulk. 



In previous sections we use the Chern-Simons approach 
to study SET phase which are non-chiral, i. e. there are no 
gapless edge excitations in the absence symmetry. Chern- 
Simons approach also applies to chiral topological phases, 
which has gapless edge modes even in the absence of sym- 
metry. These chiral phases have a nonzero chiral central 
charg#2i c_ and quantized thermal Hall effect which 
necessarily breaks time reversal symmetry. In the follow- 
ing we'll use Laughlin states as illustrative examples of 
chiral SET phases. 



C. Classifying bosonic Laughlin state at 
V = (k £ Z) with onsite Z2 symmetry 

A Laughlin state^' at filling fraction ly = 1/ni is de- 
scribed by K ~ m in effective theory ([T|) . When m =even 
it describes a bosonic topological order, while m =odd 
corresponds to a fermionic state. S uch a n effective the- 
ory also describes chiral spin liquid^^^Ell Here we start 
from the simplest case i.e. m = 2. It has 2-fold GSD 
on torus, corresponding to two different types of quasi- 
particles (or two superselection sectors): boson 1 and 
semion s. Under a unitary Z2 symmetry a semion al- 
ways transforms into a semion, hence we don't expect 
any unconventional SET phases where two inequivalent 
quasiparticles exchange under Z2 operation. Again due 
to the existence of nontrivial Z2-SPT phase of bosons in 
2+1-D, we use the following 3x3 matrix 



(31) 



in ([T]) to represent a generic = 1/2 Laughlin state with 
Z2 symmetry. The group compatibility conditions (18) 




for symmetry transformations (15) are {rj^ = 1 for uni- 
tary Z2 symmetry) 



'2 0' 
1 
.0 1 0; 



(W9)' - l3x3, I 1 I = {WBf I 1 I W9, 



'2 0' 
1 

.0 1 0; 



fl3x3+WS)(5</.S = 



/i o\ 

2 nTT, ne Z^. (32) 
^0 2 0/ 



The inequivalent solutions to the above conditions are 
W9 = 13x3 and 



They correspond to 4 different SET phases as summa- 
rized in TABLE |V] with k = 1. 

Accordingly the new quasiparticle (i.e. Z2 vortex) qg 
emerging after we gauge the unitary Z2 symmetry is 
Qg = (ji, «3, 1)^/2. Its topological spin is given by 
exp(27ri/igg), where 



27r 



2 ^ 27r 27r 



^l + 4^3 



mod 1. (34) 



Its mutual statistics with semion s = (1,0,0)^ is 



— TT mod 27r. 
2 



(35) 



Following the Chern-Simons approach to gauge the 
unitary Z2 symmetry described in AppendixjB] we obtain 
the following topological order 





(1 


^\ 

" 2 


K-i = M^K-iM, M 











hj 



and hence 








(36) 



Specifically for SET phase #4 in TABLE [Vj its 
"gauged" theory has the following quasiparticle contents 
in ([36|: 



57^ 
8 



TT, 7 = 0, 1,' 



(37) 



where (1, 0, 0)^ ~ (0, 0, 2f and (0, 1, 0)^ ~ (0, 0, 4)^. 

In a complete parallel fashion we can study generic 
"conventional" Z2-symmetry-enriched even-denominator 
Laughlin state at v = l/(2fc), fc S Z. Without loss of 
generality, a, v = l/(2fc) Laughlin state with unitary Z2 
symmetry is represented by 



K = 2fc 





(38) 



It has 12 different quasiparticles (or superselection sec- 
tors) labeled as 



97 



''q-i 



27r 2fc' 



7-0,1,. 



,2fc- 1. 



The group compatibility conditions ( 18 1 for symme 



try transformations (15) have the following inequivalent 
solutions 



f = (-y->7r,i37r)^. 



(33) 



Wf = l3x3, 



(39) 
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K ~ 2k with unitary symmetry Ga = Z2 = {g, e = g^} 



Data set in 
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Label 



On gauging symmetry: 



hqg mod 1 



mod 2-K 



[K = {2k) e 



1 

1 



,{r?'' = +l,W« = l3x3,'50n] 



#1 
(0,^,0)^ 



#2 
(0,7r,7r)^ 



#3 

(7r/2fc,7r,0)^ 



#4 

(7r/2fc,7r,7r)^ 









1/4 







/Sfc 

1 I ^ 8fc 
V 1 



_at. 



7r/2fc 




l+4fc 

aic 



7r/2fc 



Table V: Classification of "conventional' = 2^ , fe £ Z bosonic Laughlin state (or chiral spin liquid with 2A;-fold GSD on torus) 
enriched by onsite (unitary) Gs ~ Z2 symmetry. There are 4 different conventional SET phases, where under Z2 symmetry all 
quasiparticles merely obtain a (7(1) phase factor. The data set in the 2nd line completely characterizes these SET phases. Kg 
denotes the topological order, which is obtained by gauging the unitary Gs = Z2 symmetry in the Z2 spin liquid. On gauging 
the Z2 symmetry (blue entries) a new quasiparticle excitation qg is obtained. Its statistics (B3l-(B4l are also summarized in 
the table; its self statistics Og^ — 2nhqg has a one-to-one correspondence with its topological spin exp(27ri/igg). 



where ii.3 — 0, 1. The solution ii — 2 represents the same 
SET phase as ii =0, according to Corollary II in the 
criterions. Comparing with the v — 1/2 bosonic Laughlin 
state case, we can see there is a universal structure for 
all bosonic Lauglin state with K ~ 2k, k G Z. To be 
specific, for a J/ = l/2k bosonic Laughlin state, there are 
4 different Z2-SET phases as summarized in TABLE |V] 
The quasiparticles (or edge chiral bosons) transform as 



> 0+ (^,7i',«37r) 



«1.3 =0,1. 



(40) 



under "conventional" Z2 operation. After gauging the 
Z2 symmetry, one obtains a quasiparticle (the Z2 vor- 
tex) Qg = (ii, is, l)-^/2. Its topological spin is given by 
exp(27ri/i„ ) where 



1 



K 



i? + Aijk 



mod 1. (41) 



"^^ 27T 2^ 2t^ ' 2ti 
Its mutual statistics with anyon qi = (1, 0, 0)"^ is 

Oci^s = mod 27r. (42) 

Again following the Chern-Simons approach to gauge 
the unitary Z2 symmetry described in Appendix [B] we 
obtain the following topological order 



K-i = K 



M, 




and hence 




n,3 = 0,1. 



(43) 



All these 4 Abelian topological orders obtained by gaug- 
ing symmetry has | det Kg| = 8fc fold GSD on a torus. 



D. Fermionic Laughlin state a.t v — jsjqrT' ^ ^) 
with Zl symmetry is unique 

Now let's turn to the simplest fermionic Laughlin state 
with K ~ 3. It has 3-fold GSD on a torus and anyon 
excitations with statistics 6 — ±7r/3. We consider the 
following matrix in effective theory ([T]) 



K = 3 © Kt 



(44) 



where Kt generically take the form of (13 1 and (14 1. No 



tice that for a fermion system with only Z^ — {e,g = 
(—1)^^ } (fermion nu mber p arity) symmetry, there is no 
nontrivial SPT phase^^SMI 2+1-D, which hosts gap- 
less edge excitations protected by Z2 symmetry. This 
fact suggests that K = 3 is enough to describe a = 1/3 
fermionic Laughlin state with only Z^ symmetry. Now 
for such a 1 X 1 matrix K = 3, the group compatibility 



conditions (18) becomes (note that P = 2 for fermions) 



W9 = 1, 



27r 

2509 = neZ. (45) 

6 



for a unitary Z2 symmetry g. The gauge inequivalent 



solutions are 



h9 = 



with 



n G 



However notice a 



fermions in this system have gauge charge 3 in ([T]), or 
alternatively it's represented by e^"^ on the edge ([sj- 
Under the Z^ operation g each fermion obtains a 
— 1 sign, which means 3S(f>^ = tt mod 2t: and n must be 
even in (451. Therefore the quasiparticle (chiral boson) 



transform as 



a=(-i) 



<t>- 



2n+l 



n e Z. 



(46) 



According to Corollary II on smooth sewing condition 
between edges, we know different integer n e Z above 
correspond to the same Z|-SET phase. As a result when 
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only fermion number parity (Z^ symmetry) is conserved, 
the Laughlin v = 1/3 state of fermions is unique. 

It's straightforward to see that after gauging the Z; 
symmetry, we obtain an AbeHan topological order 



/ 



3 X 4 = 12. 



(47) 



In fact, the above conclusion is true for any fermionic 
Laughlin state a,t v — l/(2k + 1), k E Z with conserved 
fermion number parity. In the presence of symmetry, 
it is unique with K = 2fc + 1. After gauging the Z^ 
symmetry, we obtain an Abelian topological order = 
4(2fc + l). 



E. 



Z2 spin liquids with onsite Z2 x Z2 symmetry 




In the end we turn to a Z2 spin liquid K 



in the presence of Z2 x Z2 symmetry. The symmetry 
group Gs = Z2 X Z2 = {e,gi,g2,gig2} consists of two 
generators gi and satisfying the following algebra: 



gl=gl = {.9192? = e. 



(48) 



In an integer-spirf^ system these two generators g^ 2 can 
correspond to e.g. spin rotations along x {g^) and z (^2) 
direction by an angle of tt. Naturally the 7r-spin-rotation 
along y direction corresponds to group element 9192 ■ 

Here we'll not attempt to fully classify all Z2 x Z2- 
symmetry-enriched Z2 spin liquids. Instead, we focus on 
one nontrivial example, where spinous transforms projec- 
tively under Z2 x Z2 symmetry, in the sense that under 
27r-spin-rotation along any (a;,y, z) direction the spinon 
(or electric charge e) obtains a Berry phase —1, just like 
a half-integer spin. On the other hand, the vison (or 
magnetic charge m) transforms trivially under the spin 
rotations. Such a SET phase can be easily realized by 
e.g. Schwinger bosorP^ representation of Z2 spin liquids, 
for integer spin-S* (S* = 0, 1, 2, • • • ) 



(49) 



where a are Pauli matrices. The following constraint 
needs to be enforced for each spin 




hlh 



(50) 



to guarantee — S{S+1) for a spin-S system. Once the 
bosons form a pair superfluid (but not a superfluid) 
with {bb) 7^ (but (6) = 0), the resulting spin-S" state 
after projection into the physical Hilbert space (|50]) is a 
Z2 spin liquicP^Ml. Its spinon excitations b^ carry half- 
spin each, hence transforming projectively under 5*0(3) 
(and hence Z2 x Z2) spin rotations. By "transforming 
projectively" we simply mean that after all three sym- 



metry operations in ( 48 1 which equals identity operation 



e, all spinous obtain —1 Berry phase instead of remain- 
ing invariant (or transforming linearly). In the following 
we'll show such a Z2 x Z2 SET phase can be captured in 
the Chern-Simons approach. 

Starting from a 4 x 4 matrix Kn = ^l©(^ ^1 
to describe Z2 spin liquid, for clarity we first perform a 



GL(4,Z) gauge transformation (23 1 on Kq 



K = X'^KqX = 



/o 





-1 


i\ 








1 


1 


-1 


1 










1 





0/ 



(51) 



X = 



/l o\ 

0-10 
110 

yo 1 I J 



dctX = 1. 



We study Z2 spin liquid with Z2 x Z2 spin rotational 
symmetry in the above representation K. Notice that 



K 



K 



/ -1/2 l/2\ 

1/2 1/2 

-1/2 1/2 

V 1/2 1/2 / 



Apparently the first two components (0i.2 in the edge 
chiral boson context) can be regarded as spinous, which 
obey semionic mutual statistics with the last two com- 
ponents (</>3,4) i.e. the visons. Two spinous (visons) are 
mutually local w.r.t. each other as indicated by (|3|. 

In a Z2 X Z2 symmetry group ( 48 1 , the group 
compatibility cond itio ns for symmetry transformations 



hB^.^} in (15 



(W9i-2)^ = (W9iW92)^ = 14^4; 

(l4x4 + W9i.2)509i.2 = 27rK-ini.2; 
(l4x4 + Wfi (5092 + W92(5^9i) = 27rK-in. 

where rii, n2, n G Z^. Among all inequivalent solutions to 
these group compatibility conditions ( 48 ) , the following 
one 





/o 


1 





o^ 




/vr/2\ 




1 













7r/2 










-1 










Vo 








V 




I I 



= I4 



x4, 



A92 



= (7r/2,-V2,0,0)^ 



ni = n = (0, 0, 0, 1)^, n2 = (0, 0, -1, 0)^. 

corresponds to such a integer-spin Z2 spin liquid where 
spinous transform projectively under the Z2 x Z2 spin 
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rotations. To be specific, the quasiparticles transform as 







( 02 + 7r/2\ 






01 +7r/2 






-03 






V 04 / 






/0i+V2\ 


02 


92^ 


02 - 


03 




03 


V04y 




V 04 / 



(52) 



Indeed each spinon (0i,2) acquires -1 Berry phase after 
every 27r-spin-rotation, while visons (03,4) transform triv- 
ially. 

Notice that in such a SET phase there is no symme- 
try protected gapless edge states, i.e. generically all edge 
states are gapped in the presence of Z2 x Z2 symme- 
try. Specifically the following backscattering terms can 
be added to the edge action ([5])-(|6| without breaking 
symmetry 

^Higgs = C3 cos (203) + C4 COS (204) 

All the chiral boson modes {0i|i = 1, 2, 3, 4} on the edge 
will be gapped out by this term. 



F. Comments on symmetry protected edge states 
in SET phases 

In this section we briefly comments the symmetry pro- 
tected edge states in all SET phases discussed above. 
First of all for a chiral topological order, such as Laugh- 
lin stat^at filling fraction v = \/m,\ts edge excitations 
have net chirality |n+ — n_| ^ and hence cannot be 
destroyed even in the absence of any symmetry. These 
chiral topological orders are featured by quantized ther- 
mal Hall transport'"'"^. 

On the other hand, the edge excitations of a non-chiral 
topological order have both right and left movers and can 
be fully gapped out in the absence of any symmetry^'', 

/ 2^ 

such as in Z2 spin liquids K 



,2 0, 



and double semion 



theory K 







In the presence of global symme- 



try Gs, they might have symmetry protected gapless edge 
modes, if the edge backscattering terms are forbidden by 
symmetry. To be specific, the backscattering terms are 
typically Higgs terms ([s]). The edge states will be fully 
gapped, if and only if each chiral boson fields 0; is ei- 
ther pinned at a classical minimal by the Higgs terms or 
doesn't commute with at least one Higgs terrrPli. 

Take Z2 spin liquids with K^f*^ ^Iffil'^ ^1 for ex- 

ample, either Ai cos(20i +ai) or A2 cos(202 + (^2) could 



gap out chiral boson fields 0i,2, since [0i(a;), 02(?/)] ^ 0. 
Similarly either cos(03 -I- a^) or A4 cos(04 -I- a^) could 
gap out chiral bosons 03,4. When these terms are not 
allowed by symmetry, there could be gapless excita- 
tions on the edge protected by symmetry. The sym- 
metry protected edge modes in "conventional" (Type- 
I) Ising-symmetry-enriched Z2 spin liquids are summa- 
rized in TABLE [IT] Among the 6 different conventional 
SET phases, only #1 and #3 don't support symmetry 
protected edge modes. Phase #5 provides an interest- 
ing example. Here, both electric and magnetic parti- 
cles transform projectively under the global Z2 symme- 
try. Hence the edge perturbations Ai cos(20i -|- ai) and 
A2 cos(202 -I- a2) which attempt to condense them, are 
both disallowed, leading to a symmetry protected edge 
state. 

For "unconventional" (Type-II) Ising-symmetry- 
enriched Z2 spin liquids summarized in TABLE |III| 
there are always Ising-symmetry protected Majorana 
edge modes with central charge c = 1/2 mod 1. 

For double semion theory with K=('^ '^l©('^^l 

Vo -V V V 

on the other hand, either A^ cos(20i -I- 202 + «+) 
or cos(20i — 202 + a-) can fully gap out chiral 
bosons 01,2- Meanwhile again either ^3cos(03 + a^) or 
A4 cos(04 -I- a^) could gap out chiral bosons 03,4. When 
symmetry forbids these terms on the edge, there will be 
gapless edge excitations. The results are summarized in 
TABLE |IV| for Ising-symmetry-enriched double semion 
theories. Among the 8 different SET phases, only #1 and 
#7 don't host symmetry protected gapless edge states. 



IV. CONCLUDING REMARKS 

In summary, we have presented a general framework 
to study 2+1-D symmetry enriched topological phases 
with Abelian topological order, using the Chern-Simons 
approach. It allows us to implement generic on-site uni- 
tary (or anti-unitary) symmetry in an Abelian topolog- 
ically ordered phase in 2+1-D, to differentiate whether 
two states belong to the same SET phase or not, and to 
gauge a generic unitary symmetry and extract the resul- 
tant topological order. Based on this general formulation, 
we classify all different SET phases in a series of exam- 
ples, including Z2 spin liquids with time reversal symme- 
try (TABLE g, Z2 spin liquids with unitary Ising sym- 
metry (TABLE P and [HI|, double semion theory with 



unitary Ising symmetry (TABLE IV), bosonic Laugh- 
lin states with unitary Ising symmetry (TABLE |v]) and 
others. We also show that (odd-denominator) fermionic 
Laughlin states with only conserved fermion number par- 
ity {Z2 symmetry) is unique. Consequences of gauging 
symmetries and measurable effects, such as gapless edge 
states, are also discussed for these SET phases. 

A number of directions remain. Can the approach ap- 
plied be extended to spatial symmetries? Can we extend 



22 



this framework to symmetry enriched non-AbeHan topo- 
logical orders in 2+1-D and SET phases in 3+1-D? While 
SPT phases form an Abehan group, it is presently unclear 
if the set of SET states have additional structure. We leave 
these questions to future work. 
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Appendix A: Introduction to GL{N, Z) 

GL{N,Z) is the group of all A'' x A'' unimodular ma- 
trices. All GL{N,'L) matrices can be generated by the 
following basic transformations (i 7^ j) : 



S. 



= Sa,bi^ ~ ^ad){X ~ <5aj) + ^a,]hd - Sa,iSb,j, 

Da,b = <5a,b - '^5a,N5b,N- (Al) 



T(*'^)K will add the j-th row of matrix K to the i-th 
row of K, while S'^^'^^K will exchange the i-th and j-th 
row of K with a factor of —1 multiplied on the i-th row. 
_DK will just multiply the A^-th row of K by a factor 
of -1. KT^*'-?), K^^*'^) and KD correspond to similar 
operations to columns (instead of rows) . A subgroup of 
GL{N,'Z) with determinant -1-1 is called SL{N,Z) and 
it's generated by {T'-''^) , S'^*^^)}. 

As a simple example when N — 2, group GL{2,'Z) is 
generated by the following basic transformations: 




T 



The following results will be useful 




D 




.(A2) 




1 0' 

-n 1 , 



Appendix B: Chern-Simons approach to gauge a 
unitary symmetry 

In this section we discuss how to obtain the (intrinsic) 
topological order by gauging the unitary symmetry^'^ in 
an Abelian SET phase. We'll restrict ourselves to "con- 
ventional" SET phases, characterized by data [K, {77^ = 



-HI Wf = 1 



NxN, 



t>^\g e Gs}]- In these cases the chiral 



bosons (pi only acquire U{1) phase factors (/>/ — )■ (pj + 54>^ 
after unitary symmetry operation g € Gg- When we cou- 
ple the quasiparticles to a gauge field (with gauge group 
Gs), the following gauge flux 



r(0) 



becomes deconfined excitations in the system. Since in a 
Chern-Simons theory ([ij gauge charges are always com- 
bined with gauge fluxes by the following equation of mo- 
tion 



J 



Therefore the new excitation (emerged after gauging the 
symmetry, called Z2 vortex) carries gauge charge vector 
K(5(/)^/ (27r). We will denote a unit of this new excitation 
(the Z2 vortex) as qg. Therefore the new topological or- 
der Kg obtained by gauging symmetry must include this 
excitation in its quasiparticle contents. More precisely, 
the quasiparticle contents of topological order Kg is ex- 
panded by all the integer vectors as well as multiples of 
vectors {KS(fa/{2n)} 



^9 2Tr 



1 e Z 



N 



e Z. 



(Bl) 



where {g} are the generators of symmetry group Gs that 
are gauged. And we can identify the new matrix Kg 
which contains all these quasiparticles in its spectra. 

In the following we work on one example to demon- 
strate this gauging procedure. We consider SET 

and symmetry transformation W^' — I4X4, Scf)^ = 

7r(ii/2, «2/2, 1, 44)"^ with ii.2,4 = 0, 1. According to (Bl I 



phase #6 in TABLE 



II 



with K 



we know here a generic quasiparticle is labeled by gauge 
charge vector 



1' = Ml; 



M 



fi2/2 1 0\ 

ii/2 1 

^4/2 1 

Vl/2 0/ 



lez^ 



Since this new Abelian topological order is determined by 
the statistics of its quasiparticles, we immediately obtain 



(i;)-K-4;, = ifK-ii2, 1^ 



ML 



and therefore 



K-i = M^K- 



/iii2±2U ,^/4 ,2/4 i/2\ 



iM 



K„ 



h/4 
V 1/2 

fo 

2 

2 

\2 -ii 



1/2 
1/2 


-42 

~2i^) 





/ 



(B2) 
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Clearly from (12|-(|3| we know the self statistics of new 
quasiparticle qg is 



501 _ _ M2+2i4 
27r ~ " 4 



(B3) 



and its mutual statistics with original quasiparticles 
e, m, / of Z2 spin liquid are 



J = 



2 ' 



'Ig, 



2 ' 



q : n _ (n+i2)7r 



(B4) 



Topological spirpSl exp(27rift,q). the Berry phase ob- 
tained by adiabatically rotating a quasiparticle q by 27r, 
is an important character of a 2+1-D topological or- 
der. In Abelian topological orders, the topological spin 
exp(27rifeq ^) h as a one-to-one correspondence to the self- 
statistics (|B3| of a quasiparticle in unit of 27r: 



"Ig 



"Qg 

27r 



h9. 



1 



2^ 27r 



iii2 + 214 



27r 



(B5) 



For the "unconventional" SET phases, e.g. in our case 
with Gs = Z2, the Z2 symmetry would exchange quasi- 
particles that belong to different superselection sectors 
in Abelian topological order K. Gauging this kind of 
Z2 symmetry will in general lead to U{1)^ x Z2 Chern- 
Simons theory which describes non-Abelian topolog- 
ical orders in relation to Z2 orbifold conformal field 
theorjE2l£Hl. 



Appendix C: Classifying double semion theory with 
onsite Z2 symmetry 



Double semion theorji^^lMI jg a "twisted" Z2 gauge the- 
ory in 2+1-D. with Abelian topological order described 
^2 0^ 

bosonic Z2-SPT phase in 2+1-D, again here we believe a 
4x4 matrix 



by K = 



Due to the presence of nontrivial 



K 




/2 






o\ 
-2 
1 
1 0/ 



(CI) 



can capture all the different Z2-symmetry-enriched dou- 
ble semion theory. Such a theory has the following quasi- 
particle contents in its spectra 



(C2) 
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v) 
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vo/ 



where s and s represents semion and anti-semion re- 
spectively, and b is the bound state of a semion and an 
anti-semion. b has bosonic self statistics ([2| but mutual 
semion(anti-semion) statistics with s(s). Here {l,s,s, 5} 
represents the 4 superselection sector of double semion 
theory. Any two quasiparticles differing by a local exci- 
tation ~ belong to the same superselection sector. 

Now let's consider the implementation of unitary Gg — 
Z2 symmetry on double semion theory. We have group 
compatibility condition ( 18 1 for symmetry transforma- 



tion (15 1 on quasiparticles: 
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W9, 



n, n G 



We consideilSSI the solutions to ([CS]) with = 14^4. 
Due to Cri terion I, the gauge inequiva lent solutions of 
6(j)B to ( C3 ) are summarized in TABLE 



IV 



Following Appendix^ we briefly discuss consequences 
of gauging the unitary Z2 symmetry in the double semion 
theory. Since the symmetry transformation is a C/(l) 
phase shift Scfi^ — 7r(zi/2, 22/2, 1, m)"^ as shown in TA- 
BLE |IV] the quasiparticle content in the new topological 
order obtained by gauging Z2 symmetry is expanded by 
gauge charge vector: 



1' = Ml; M = 



/ ii/2 1 0\ 

-^2/2 1 

^4/2 1 

V 1/2 0/ 



le 



and therefore 



f^:!^^ *i/4 Z2/4 l/2\ 
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Take #7 as an examples with 5$b = 7r(l/2, 1/2, 1, 0)"^ 
(or ii = 12 = 1,^4 — 0), the Abelian topological order 
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obtained by gauging Z2 symmetry is 

/O 4 0\ 

4 



K„ 



K„X = 



X = 
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Vo 1 oy 

/ 1 o\ 

110 
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\-2 2 11/ 




e GL(4,Z). 



Notice that 




Again from ([2])-([3]) we can obtain the self statistics of new 
quasiparticle qg as 



^98 



2^—^2+424 



(C5) 



and its mutual statistics with original quasiparticles 
e,m, f of Z2 spin liquid are 



9qgM = 



2 ' 



_ 12TT 
~ 2 ' 

(ii+i2)7r 



^..rs (C6) 

The topological spin of this new quasiparticle is given by 
exp(27ri/ig_^) where 



"la 



1 .5^sT 

2 ' 27r ^ 27r 



*1 ~ *2 + 4«4 

16 



(C7) 



Unlike others, for SET phases 7^4, 7^6, #8 it's not easy 
to find a GL{A, Z) transformation (111 on Kp matrix ( C4 ) 



to reduce it to a simpler form. e.g. one can only show for 
SET phase #4 



X^KgX 
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and for SET phase #6 
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However, a one-to-one correspondence between the quasi- 
particle contents of two seemingly different Kg matrices 
can be established. For example, there are 16 differ- 
ent superselection sectors (or 16 quasiparticle types) for 
Abelian topological order Kg in ( C4 1 , obtained by gaug- 
ing Z2 symmetry in ^4 SET phase (ii = 1 = 24,12 = 0): 

/71 + 472\ 





7i + 72 




in (C4l ^ 




71 =0,1,--- ,7, 72 = 0,1 
For #6 SET phase (^2 = 1 = 14, n = 0): 
/ 471 -I- 372 \ 



71 + 72 





in (C4l ^ 




72 =0,1,--- ,7, 71=0,1. 

The Abelian topological order obtained by gauging Z2 
symmetry in #8 SET phase (12 = ii ~ ii ^ I) , however, 
doesn't correspond to any of the usual Z2 x Z2 (or Z4) 
gauge theory. Instead it is characterized by the following 
matrix 
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X^KgX 



(0 


2 





o\ 


2 


2 


1 








1 


2 


2 


\o 





2 


oy 



(C8) 



/O 



X 



Vo 



detX = 1. 



1 l\ 

-10 


1 oy 

And it has 16 different types of quasiparticles: 
^\ 

71,2 = 0,1,2,3. 



7 ^ 



72 



voy 



in (C4l 



Among them four has bosonic self statistics 
mod 2tt), six with semionic statistics {6 = ^ 
and the other six with anti-semionic statistics 



(6* = 
mod 2tt) 



mod 2 tt). I n the above basis, the 16 x 16 modular 5- 
matrijPsESl of this Abelian topological order is given by 



■ exp 



2 

Y(2 5l7a7; 

a=l 



+ 7i72 + 727i) • (C9) 



Appendix D: Vertex algebra approach to gauge a 
unitary symmetry 

The Chern-Simons approach to gauge a unitary sym- 
metry, introduced in Appendix |B] applies to all cases 
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where we obtain an Abelian topological order after gaug- 
ing the symmetry. Thus for any "conventional" SET 
phases we can gauge its unitary symmetry and obtain 
an Abelian topological order in the Chern-Simons ap- 
proach. For "unconventional" SET phases, such as those 
summarized in TABLE III gauging a unitary {e.g. Z2) 
symmetry will result in non-Abelian topological orders. 
In the case Gg = Zi as discussed in this work, these non- 
Abelian topological orders are described by [/(l)^ x Z2 
Chern-Simons theory"^. In these "unconventional" cases 
the Chern-Simons approach introduced previously is not 
enough. In order to obtain the full structure (such as 
topological spin exp(27ri/i) of quasiparticles and modular 
S matrix associated with quasiparticle statistics) of these 
non-Abelian topological orders, here we introduce a ver- 
tex algebra approach to gauge the unitary symmetry. It 
applies to both the "conventional" and "unconventional" 
SET phases and in the following we'll demonstrate its 
power by two examples: "conventional" and "unconven- 
tional" Z2-symmetry-enriched Zi spin liquids. 



The vertex algebra formalism, and application 
to "conventional" SET phases 



The vertex algebra approadP^EsEI] ig based on the 
close connection'*-^* between the bulk topological or- 
der (described by 2+1-D topological field theory) and 
its boundary excitations (described by 1+1-D conformal 
field theory) in two spatial dimensions. Let's take Z^ 
spin liquid ( [24| for an example. The edge effective theory 
([5]) contains two branches of chiral bosons {(/)i,2}, which 
could be reformulated by ac = 1 C/(l) x J7(l) Gaussian 
model with a holomorphic and anti-holomorphic part: 

Lp(x + ir) = ip{z) = 01 (x, t) -f (/)2(x, t), 
ip{x - ir) = ip{z) = (t)i{x,t) - (f)2{x,t). 

The Gaussian model has Lagrangian density Coaussian = 
^dip{z)dip{z) = g^|V(^p, yielding the following correla- 
tion function 



{ip{z)(p{w)) = — ln(z — w). 



(Dl) 



and {(p{z)(p{w)) ~ — ln(z — u)). The free boson field (p has 
compactification radius R = 2 for Z2 spin liquid (24) so 
that periodicity ip 

f N 



in (5 



\N 0^ 
scalar boson 



2ttR holds. In general for K = 
the associated compactification radius of 



jrViz) = (t)l{x,t) + (f>2{x,t), 



(p{z) = 4>l{x,t) - (f)2ix,t). 



(D2) 



is i? = \/2N. The allowed physical excitations must be 
compatible with 27ri? periodicity of bosons and they are^ 



for holomorphic part (and similarly Vk for anti- 
holomorphic part). These 27V vertex operators are pri- 
mary fields of the holomorphic U{1) conformal field the- 
ory (CFT) and they form different representations of the 
conformal algebra. From (Dl I one can see they have the 
following (radial-ordered) operator product expansioEp^ 
(OPE): 

giav(z)gi/3yM u;)"/3e'("+'3)'^('") + . . . (D4) 

for a + f3 7^ 0. There is an energy- momentum tensor T — 
— ^{dip)"^ which generates the conformal transformation 
of the vertex algebra, so that any primary field P{z) has 
the following OPE with energy-momentum tensor 

T{z)P{w) = , „ P(w) + —^dP{w) + ■■■ (D5) 
(z — wy z ~ w 

where hp is the scaling dimension of primary field P. 
Apparently the vertex operator exp [1095(2)] has scal- 
ing dimension — \cP' ■ Another primary field is the 
current operator j{z) = id(p{z) which have scaling di- 

aj{w) 



mension hj = 1. And we have 



{z-w)°'^ (z-w)""-! 



These OPEs imply the following fusion rules of primary 
fields 

e'"^ X e'°"P = 1 + j, (a^ 0) 

Similar results hold for anti-holomorphic (p{z) part, only 
that all scaling dimensions changes sign for their anti- 
holomorphic counterparts. 

A natural question is among all these primary fields, 
which ones appear in the physical edge spectra of the 
topologically ordered phase? There are a few physical 
principles to follow. First of all every physical edge exci- 
tation (e'^i'*"^', k € Z) must be a primary field. Sec- 
ondly, there are electron operators (or the microscopic 
local degrees of freedom e' j i.iKi.j<l>j ^ g which is 
local with respect to all other edge excitations. In the 
context of vertex algebra, two operators A and B are 
local w.r.t. each other if and only if in OPE 

A{z)Biw) = ^^P^ + 0[{z-wy-'^---), (D6) 

where C is also a primary field and f^g is a structure 
constant. For example before we gauge the Z2 symmetry, 

^0 

0^ 

operator is 



for Abelian topological order K 



the electron 



exp 



Vk{z) 



0,1,- •• ,27V- 1. (D3) 



h.2 e z. 



(D7) 
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It's straightforward to check that all allowed quasiparti- 
cles (local w.r.t. the above electron operator) have the 
following form 



3i(ilI^l+i2</>2) 



exp 



ih+l2)v{z) + {h-l2)v{z) 
V2N 



h,2 e Z. 



Lastly, any two primary fields differing by an electron 
operator are regarded as the same (or belong to the 
same superselection sector). 



Now let's go back to Z2 spin liquids with K = 



which have 4 branches of chiral bosons {<f)i,l < 





2. Application to "unconventional" SET phases 

For a "unconventional" SET phase, e.g. where two in- 
equivalent quasiparticles (e and m in Z2 spin liquid) are 
exchanged under Z2 symmetry operation as summarized 
in TABLE |III| a non-Abelian topological order is ob- 
tained by gauging the Z2 symmetry. Here we apply the 
vertex algebra approach to extract the full structure of 
these non-Abelian topological orders. 

First let's review some known results, discussed in de- 
tail in Ref. 155)721 When the "unconventional" Z2 sym- 

^ ^ 6^ = 0} is gauged for Abelian 



metry {W^ 



1 0, 



We can introduce free bosons ipi{z), ipi{z) for topological order K 



the resultant topologi- 



chiral bosons 0i_2 as in (D2 ) with N — 2, an d free bosons 
ip2{z) , (p2{z) for chiral bosons ^3^4 as in (D2 1 with N 
In other words we have 



1. 



fl/2 
1/2 

V 



1/2 
-1/2 



I/V2 



I/V2 
-1/a/2 



Before gauging the unitary Z2 symmetry, the four super- 
selection sectors (or 4 types of different quasiparticles) 
correspond to 



3i{z) - ji(z) - e 



2ii^i(z) 



■ j2{z) ~ 



n{z) 



%/2 



m ' 



yi(^)-'y^l(g) 



Now after gauging the "conventional" Z2 symmetries in 
TABLE |IT] as discussed in Appendix [B] a new type of 
quasiparticles qg becomes deconfined excitations: 



' N 
0, 

cal order is described by U{1) x U(l) x Z2 Chern-Simons 
theory (coined "twisted" Z^ gauge theory in Ref. 17^ . 
which has GSD= {NS/2)[NS -|- 1 + {2^9 _ l)(Ars-i _^ ^-^^ 
on a genus-5 Riemann surface. It contains 2N different 
quasiparticles with quantum dimension d = 1, another 
27V quasiparticles with d = y/N and N{N — l)/2 quasi- 
particles with d = 2. Under the unconventional Z2 sym- 
metry operation two superselection sectors e m ex- 
changes and so does chiral bosons 0i <^ 4>2- Therefore in 
the context of vertex algebra (D2) the anti-holomorphic 
free boson (p ^ —(p under Z2 symmetry operation! After 
this "unconventional" Z2 symmetry is gauged for Z2 spin 
liquids {N = 2), we obtain an non-Abelian topological 
^211^2(2) ^ gv^i 1^2(2) order whose quasiparticle content has an antiholomor- 
phic part (from (f) given by Z2 orbifold CFT with 
compactification radius R = 2. It has been shown that 
(I)8)Z2 orbifold CFT is equivalent to Isingxising (or Ising^) 
CFT^. In each Ising CFT there are 3 different quasipar- 
ticles: vacuum (or boson) 1, fermion ip and the "disorder" 
field^* a with the following fusions rules: 



,i \ viiz)±Vi{z)\ 



jiVl(z) 



1:1. 



(»i+»2)yi+('2-'i)'Pi+y2(i+'4)y2+V2(u-i)y2 



exp 



where ^1,2,4 = 0, 1 in 6(j)^ . Notice that when such a Z2 
vortex Qg is deconfined, we have to modify the previous 
definition of electron operators 1 in (D8I. The new elec- 
tron operator is defined as anything that is local w.r.t. 
quasiparticles {e,m, f,qg}. With this new definition for 
electron operators, we can track down all the inequiva- 
lent quasiparticles (superselection sectors) and obtain the 
full structure of the topological order obtained by gaug- 
ing Z2 symmetry. One can easily check this approach 
indeed reproduces TABLE |TT] consistent with the result 
of Chern-Simons approach. 

In the vertex algebra context, the scaling dimen- 
sion h oi a quasiparticle determines its topological spin 
exp(27ri/i), the Berry phase obtained by self-rotating a 
quasiparticle adiabatically by 27r. On the other hand, 
the mutual statistics of quas ipart icle A and B is given 
by 0A,B = -27raA,B in OPE (|d6|. 



^ X tp ~ 1, tp X a = (7, C7 X a = 1 + 'ip. (D9) 

Both 1 and -0 have quantum dimension 1 while disorder 
operator a has quantum dimension \/2- Their scaling 
dimensions are 0, 1/2 and 1/16. Therefore the Z2 
orbifold CFT, equivalent to the direct product of two 
copies of Ising CFTs, contains 9 = 3x3 inequivalent 
quasiparticles (superselection sectors). The quasiparticle 
contents of the Z2 orbifold CFT are summarized in the 
first 3 columns of TABLE IVll 



Now let's get back to our cases of Z2 spin liquids with 
unconventional on-site Z2 symmetry. There are 4 such 
SET phases as summarized in TABLE |TTT] After gauging 
the unitary Z2 symmetry, they all lead to non-Abelian 
topological orders with 9 inequivalent quasiparticles (su- 
perselection sectors). In the vertex algebra context, they 
all share the same antiholomorphic (ipi) part which gives 
rise to the non-Abelian quasiparticles. However, their 
different holomorphic parts discriminates these 4 SET 
phases. A key issue in determining the quasiparticle con- 
tents is: which quasiparticles are identical (or belong to 
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Relation to Kitaev's 16-fold wajBSl 




{u = 1) (g> {ly = 15) 




(// = 5) ® (i^ = 11) 
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(iy = 3)(g){u = 13) 



Table VI: Quasiparticle (q.p.) contents of non-Abelian topological orders obtained by gauging the Z2 symmetry in "unconven- 
tional" SET phases as summarized in TABLE |III[ They are related to Z2 orbifold CFT compactified at radius R — 2. The 
fusion rules of non-Abelian quasiparticles have a one-to-one correspondence to the two copies of Ising CFT (i.e. Ising^ theory). 
Each quasiparticle qa correspond to a vertex operator (a primary field) in the vertex algebra (which are CFTs) defined through 
operator product expansion (OPE), and its (conformal) scaling dimension ha ( mod 1) physically relates to the topological 
spin exp{2niha) of the quasiparticle. The modular 5 matrix of such non-Abelian topological orders is also determined by the 
OPEs between vertex operators. Allowed quasiparticles must be local w.r.t. any electron operators ~ 1. Any two quasiparticles 
differing by an electron operator are considered as the same. The scaling dimensions in the Ising^ CFT (or Z2 orbifold model) 
are hi — 0, hj — 1, hfi — hf2 = 1/2, hv\ = 1/8, h^i = ha.2 — 1/16 and h^i = h^2 = 9/16. We label these 9 different 
quasiparticles (or superselection sectors) as g^, < a < 8, which is shown in the (a-|-2)-th row of this table. All 4 non-Abelian 
topological orders (7^1 — #4) have 9-fold GSD on a torus, corresponding to 9 different superselection sectors. 



the same superselection sector), after the symmetry is 
gauged? 

In the vertex algebra context, once we fix the elec- 
tron operator 1 ~? (or the trivial sector) which is local 
w.r.t. all quasiparticles, the full structure of inequivalent 
quasiparticles is determined. So the above issue becomes 
the following question: how to determine the electron 
operators in the vertex algebra, once we gauge the uni- 
tary symmetry? The answer lies in the following physical 
principle: 

// in the original SET phase, two quasiparticles belong 
to the same superselection sector (i.e. they are equivalent) 
and transform in the same way under a unitary symme- 
try, then they belong to the same superselection sector 
after the unitary symmetry is gauged. 

To be specific, if two quasiparticles qA and qB belong 
the same superselection sector and transform in the same 
way under unitary symmetry, then after gauging the sym- 
metry, quasiparticle qAQs ~ ^ (Ib anti-particle of 
qs) belong to the trivial sector. For instance, in SET 



phase #1 in TABLE III and | VI| the following two quasi 



particles belong the the trivial sector and are both odd 
under Z2 symmetry g: 



Jl 



• VJ2 + 'P2 

e ^ 



and they are both their own anti-particles. Besides the 
following two fermions also belong the same superselec- 
tion sector and are both even under Z2 symmetry: 



P 



02) = C0s((/3i) 



gi(0i-H02) _ g 



Both of them are also their own anti-particles. Therefore 
we have the following definitions of electron operators (or 



trivial sector) as shown in TABLE VI 



1 jie 



iipi 



This enables us to obtain all the 9 inequivalent quasi- 
particles (superselection sectors) as summarized in TA- 
BLE I VI I for the non-Abelian topological order acquired 
by gauging Z2 symmetry in these SET phases. 

As discussed earlier, in the vertex algebra approach, 
the mutual statistics of two quasiparticles A and B is 
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given in their OPE (D6l by Sa.b = gxpCiOa.b) = 
exp(—2'KiaA,B)- If quasiparticles A and B leads to 
more than one fusion channels, the corresponding entry 
Sa,b — vanishes in the modular S matrix. Besides, 
scaling dimensions {hq} of quasiparticles {q} determine 
their topological spins Ta,b — <5a,b exp(27ri/i/i), which 
corresponds to the modular T matrix. So we can extract 
all the topological properties of the non-Abelian topo- 
logical orders, obtained by gauging Z2 symmetry in SET 
phases. 

The modular S matrix in the basis Qa {0 < a < 
8, see TABLE IVip of the 9 different quasiparticles (su- 
perselection sectors) is S^i — 
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for gauged "unconventional" SET phase #1. Mean- 
while the T matrix is a diagonal unitary matrix Ta,h = 
(5a,b exp(27ri/ia), where ha gives the topological spin 
exp(27ri/ia) of quasiparticle qa shown in TABLE VI To 
be specific we have 
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-i7r/8 J 



For SET phase #2, after gauging the unitary Z2 sym- 
metry we have its modular S matrix as S^2 = 
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and its T matrix as 
/I 



T#2 
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,i5ir/8 
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,i57r/8 
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Clearly after gauging the unconventional Z2 symmetry, 
SET phases #1 and #3 lead to the same non-Abelian 
topological order. They share the same S and T matri- 
ces, differing by a relabel of quasiparticles in TABLE [VTl 
For example quasiparticle (75 in ^1 phase corresponds to 
quasiparticle in #3 phase. Similarly SET phases #2 
and #4 lead to the same non-Abelian topological order, 
by gauging the unconventional Z2 symmetry. 

It's easy to verify that they satisfy the following con- 
sistency conditionals! for modular transformations: 



1. 



{srf 

where the C/(l) phase factor 9 is defined as 



^g2xic_/8_ 



(DIO) 



(Dll) 



da and ha corresponds to the quantum dimension and 
topological spin exp(27ri/iQ) of quasiparticle qa respec- 
tively. c_ is the chiral central charge of the edge exci- 
tations of the topological ordered phase. All the non- 
Abelian topological orders in TABLE VI have c_ = 
and hence = 1. In fact for all the 4 non-Abelian topo- 
logical orders summarized in TABLE VI their 
modular S and T matrices satisfy = [ST)^ = Igxo- 

Starting from a Z2 gauge theory {Z2 spin liquid or dou- 
ble semion theory) with unitary Z2 symmetry, once the 
symmetry is ga uged, a resultant Z2 x Z2 gauge theory 
is expectecP^EIl xhe above non-Abelian topological or- 
ders can be regarded as "unconventional" Z2 x Z2 gauge 



theories, related to Kitaev's 16-fold way classification^Sl 
of Z2 gauge theories in 2+1-D. In particular, they are 
associated with Z2 gauge theories where fermions hav- 
ing an odd Chern number {v =odd) couple to Z2 gauge 
fields. Notice that before gauging the symmetry, all 4 
SET phases have non-chiral edge excitations with chiral 
central charge c_ =0. As a result, we expect that after 
gauging the Z2 symmetry their edge states remain non- 
chiral and should be gapped due to backscattering in a 
generic situation. Indeed in all the "gauged" non-Abelian 
topological orders in TABLE |VI| a Z2 gauge theory with 
fermion Chern number v is always accompanied by its 
time-reversal counterpart v = \Q — u mod 16 through a 
direct product. 
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To be specific, in Ref. I^HlKitaev introduced a 16-fold 
way classification of 2+1-D Z2 gauge theories, describing 
fermions coupled to a Z-i gauge field. When the Chern 
number v of fermions changes by 16, one ends up with 
the same Z2 gauge theory. Specifically when v =odd, as- 
sociated Z2 gauge theory contains 3 inequivalent quasi- 
particles: vaccum (or boson) 1, fermion i\} (e in Kitaev's 
notatiorP Sl) a nd vortex a. Their fusions rules are the 
same as (D9l, i.e. those in Ising anyon theory^^. Their 



quantum dimensions are 



di = = 1, da = V2. 



Therefore when 1^ = 1 this corresponds to the Ising anyon 
theory. When a direct product of a Z2 gauge theory with 
Chern number i/ (we denote this Z2 gauge theory by i/) 
and its time reversal counterpart i? = 16 — v is made, 
one can combine the fermion ip in i' and the vortex a 
in D to form a new vortex operator, which have scaling 
dimension \ — = Therefore one can clearly see 

the following two seemingly different direct products 



J/ (g) (16 - z/) ~ (8 - i^) (g) (8 v). 



(D12) 



The topological spin exp(27ri/i) of these quasiparticles 
are given by 



hi =0, = 7:, ha = 



16' 



lead to the same topological order. As a result SET 



phases #1 and #3 (#2 and #4) in TABLE |III| lead to 
the same non-Abelian topological order, by gauging the 
unitary Z2 symmetry. 
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Although the formulation of (7(1)^ Chern-Simons theory 
in |T]) seems to suggest existence of A'^ conserved (7(1)- 
currents, they can actually be explicitly broken by e.g. in- 
troducing Higgs termJ^. These Higgs terms merely con- 
dense local bosonic excitations (instead of nonlocal anyonic 
excitations) and hence don't change the superselection sec- 
tors or any topological properties of the phase. 
In fact there are other solutions to conditions (271, such 

as W = -l4x4, W = -a, e l2x2, W = ±l2x2 e 

and — ±cjx © u^, where (Jx,y,z are the three Pauli 
matrices. Whether they may lead to new SET phases are 
is clear to us at this moment, and we don't consider them 
in this work. 

In a half-integer spin system, on the other hand, spin 
rotations by an angle of 2-n will lead to a Berry phase 
— 1. Therefore the group structure generated by 7r-spin- 
rota tions along x and z directions is not Z2 x Z2 as in 
(48|. 

We believe there is no "unconventional" implementation of 
onsite Z2 symmetry in the double semion theory. This is 
because here the three types of anyons s, s, 6 have differ- 
ent self statistics. As a result exchanging any two of them 
shouldn't be a symmetry of the system. 



